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1. Introduction 

Can quantum mechanics, which determines behavior at the atomic and sub- 
atomic scales, be manifest at a macroscopic scale? This question, which was posed 
when the foundations of quantum theory were first laid, has fascinated physicists for 
more than seventy years. 

The phenomena of superconductivity and superfluidity in helium, are quantum 
manifestations on a macroscopic scale. In both cases there is a macroscopic non- 
dissipative current of particles. More recently, quantum manifestations at scales well 
above the atomic scale were observed - for example, quantum tunneling of the phase 
in a Josephson junction, permanent currents in small conductor rings and more 

recently, Bose condensates. These systems, whose sizes vary from 10 to 10^ nm, are 
relatively complex; nevertheless, their properties can be described using a small 
number of degrees of freedom defined as a "macroscopic order parameter". 

In magnetism, since the discovery of superparamagnetism by Neel, it has been 
known that a ferromagnetic or ferrimagnetic particle of a few nm in size can also be 
described with a small number of degrees of freedom, those of the magnetic moment 
of the single domain particle, which behaves as a small magnet (the exchange energy 
dominates by orienting all the moments in one direction). 



The search for quantum effects at the "macroscopic scale" in magnetism started 
in the early seventies after it was shown that single crystals of rare-earth intermetallics 
(Dy 3 Al 2 , SmCo3.5Cui.5) exhibit fast magnetic relaxation in the Kelvin range. This 
phenomenon was interpreted in terms of magnetization reversal by quantum tunneling, 
below a certain crossover temperature. The magnetization reversal of the bulk crystals 
being the sum of the elementary reversals of single domain blocks (nucleation of the 
so-called Barkhausen jumps), this type of study is equivalent to the study of single 
nanoparticles, but with many complications (size, energy barrier and switching field 
distributions, effects of domain walls, various dissipation effects...). 

More recently there have been developments in various disciplines that have 
led to great progress in the obtention of different types of nanoparticles. In material 
science, magnetic materials have been produced as isolated aggregates, as deposits of 
aggregates, as carbon nanotubes and nanocages filled with magnetic material, as 
electrodeposits of magnetic material in nanoporous polycarbonate membranes, as well 
as dispersals in polymers. Molecular chemistry has produced molecules with giant 
spins and colloidal chemistry has used micelles as microreactors to make all sorts of 
new magnetic nanoparticles. Naturally occurring biological systems have given us 
ferritin and biochemistry has provided us with their artificial analogues. 

Among all these systems, an exciting type of material has emerged for the study 
of macroscopic tunneling in magnetism - namely, molecular crystals with identical 
magnetic molecules. Despite their relatively large size and the dipolar interactions 
between their magnetic moments, these molecules clearly exhibit quantum tunneling of 
magnetization. The focus of this article is on two of these materials (i) the spin cluster 
system so-called "Mni2-ac", which is a Manganese acetate. This system was the first to 
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exhibit what is referred to as "resonant tunneling of magnetization", (ii) Another 
analogous system is the so-called "Feg", in which the same phenomenon was observed 
afterwards, but at lower temperature, which allows easier experimental studies. 

2. Mnn-acetate. 

2.1. Experimental results. 

Over the past years a lot of experimental and theoretical works were performed on 
molecules of [Mni20i2(CH3COO)i6(H 2 0)4]. This molecule has a tetragonal symmetry 
[1] and contains a cluster of twelve Mn ions divided into two shells (four Mn 4+ ions 
from inner shell with spin S=3/2, surrounded by eight Mn 3+ ions from the outer shell 
with spin S=2) with strong antiferromagnetic couplings (frustrated triangles, see Figure 
1). They form a collective ground state spin S=10 with magnetic moment M=gS|i B 
~20(l B (g~2 is the Lande factor) [2]. These molecules are chemically identical and 
form a crystal with an average distance between Mni 2 molecules of the order of 15 A 
[1]. Intermolecular exchange interactions are negligible and dipolar interactions are 
about 0.01-^0.02 K. This is much smaller than the anisotropy barrier U of each 
molecule (which is about 61 -h65 K [6(a),7, 9,10]). 

Ac-susceptibility measurements of Sessoli et al. [2] and magnetization 
experiments of Paulsen et al. [6(a),7] show a superparamagnetic behavior with a 
relaxation time that obeys the Arrhenius law T=T exp(U /k B T) (see, L. Neel in [3]) with 
x o =2xl0" sec at high temperatures (T>2.5 K) and a blocking temperature T B close to 3 
K (TB~k B Uo/ln(t/x )~3.3 K for t~l h). Above the blocking temperature, this 
superparamagnetic behavior is characterized by a Curie- Weiss law with very small 
positive paramagnetic temperature 0-70 mK. This indicates the existence of weak 
dipole-dipole interactions between molecules. As the temperature decreases from ~3 K, 

3 



the magnetization evolves from a relatively fast roughly exponential relaxation (about 
10 3 sec) to very slow non-exponential relaxation ( about 10 7 sec at T-2.1 K ; a non- 
exponential behavior discovered later, will be discussed below) [6(a),7,9] (experiments 
[8] even give logarithmic relaxation below 1 K, but this not confirmed by more recent 
results [*](I. Chiorescu, R. Giraud, A. Canneschi, L. Jansen, and B. Barbara, submitted 
for publication). It was also found that the relaxation time x (H) exhibits a deep 
minimum in zero field at T<2 K (see, for example, [7,9-11]) whereas above this 
temperatures x (H) presents a maximum at 0.2 T (see Figure 2) observed, for the first 
time, by Paulsen et al. [7] (see, also, review paper of Barbara et al. [9]). This picture 
(strongly supported by dips observed in ac-susceptibility measurements by Novak et 
al.[10]) was interpreted as Quantum Tunneling of Magnetization (QTM) of the 
collective spin S=10 with a crossover temperature T c ~2 K, due to the resonant energy 
level crossing in two- well potential in longitudinal magnetic field [6(a),7,9-ll,4]. 

The experiments reported in [6a, 9] were thoroughly repeated and confirmed in 
[7]. More detailed measurements performed later [4,11], showed steps in isothermal 
hysteresis loop (Friedman, et al. in [4] and Thomas, et al in [11]) when the field is 
increased in the direction opposite to the magnetization at T<T B (Figure 3, see, also 
[6(b)]). In the flat regions of the hysteresis loop the relaxation times were found much 
longer than the experimental time window (-600 sec in [11]) while in the steep regions 
magnetization relaxes much faster and the relaxation times can be of the order (or even 
less) of the experimental time window. The plot 3 M z /3 H z versus longitudinal filed H z 
(Figure 4) gives a series of peaks (with Lorentzian shape). The maxima of these peaks 
allows to define the values of magnetic field H n at which the magnetization steps occur : 
the fields H n such as H n ~0.44n T (n=0, 1, 2, ...). The relaxation measurements of 
Thomas et al. [11] show that the relaxation time oscillates with the magnetic field 
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(Figure 5) and has deep minima (resonances) at the same values of field where the steps 
are observed in the hysteresis loop. These new experimental results clarify the problem 
of the maximum observed before at H=0.2 T [7,9,10]. As suggested by Barbara et al [9], 
it is effectively the first maximum of x(H) curve. In addition, the blocking temperature 
T B also exhibits strong minima approximately at the fields H n (extracted from the 
temperature behavior of magnetization at different values of field in [4]). 

All these results, interpreted as a strong evidence of QTM, were obtained at 
T<T B = 3.3 K (in a longitudinal field). At the higher temperatures the magnetization 
relaxes too fast for quasi-static measurements and ac-susceptibility 
Oc(co)=%'(u))- i%"(co)) measurements are necessary. The relaxation times can then be 
determined, either from the position of the maxima of the imaginary susceptibility 
(%(u))" is maximum at cox=l) or from the relation x=%"(to) /co(%'(co) - %'(°°) ), where 
(Ois the frequency of the ac-field. The higher temperature relaxation times also 
oscillates with minima at nearly the same field as in the low temperatures case, but the 
amplitude of oscillations decrease. At e.g. 10 K this amplitude is 25 times smaller than 
in the low-temperature case [12] (Figure 6). These measurements show that, together 
with the regular decrease of the relaxation time with the applied field (usual field 
dependence of the barrier), there is even at temperatures well above the quasi-static 
blocking temperature some fraction molecules with tunneling channels, as in the low- 
temperature regime. This means that at high temperature, the mechanism of relaxation is 
intermediate between the quantum regime (where the ratio of the relaxation times at 
resonance and out of resonance are much smaller than 1) and the classical regimes 
(where this ratio goes to 1). 

It was also experimentally found in [4,11] that the transition rate decreases rapidly 
with temperature. Therefore, if the QTM mechanism is relevant for all these 
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experiments, then it can be understood only by Thermally Assisted QTM (introduced by 
Novak and Sessoli in [10] and Barbara et al. in [9]), where tunneling occurs from 
excited levels. As the temperature decreases, the higher levels (close to the top of the 
barrier) become less and less populated and the tunneling takes place between the lower 
levels with smaller probability, explaining the fact of the extremely long relaxation time 
observed below 1 K in [6(a),7,9]. 

However in the presence of increasing longitudinal or transverse magnetic field, 
the relaxation is faster and can be easily measured in the main bulk phase of Mni 2 (see 
very recent experiments [85]). As an example we show in Figure 7(a) some hysteresis 
loops obtained from torque experiments in fields up to 6 T. As in Figure 3, where 
measurements were done in lower fields and higher temperatures, the loops depend on 
temperature, but here this is true only above 0.8 K. Below this temperature the loops are 
independent of temperature, showing that tunneling takes place from the ground state 
S=10. 

Similar results were obtained in a transverse magnetic field of about 4 T. As an 
example we give in the insert of Figure 7(b) the temperature dependence of the 
relaxation times measured at short time-scale. They are clearly independent of 
temperature below 0.8 K (this temperature is nearly the same as in longitudinal field by 
accident). This result shows that tunneling occurs here between the ground states S=10 
and S=-10 of the two symmetrical wells. The main part of Figure 7(b) gives also an 
example of relaxation curves measured in the plateau. At short time-scale the relaxation 
follows a square root law, but at long time-scale it is exponential (we should note that, 
although the two third of the relaxation curves could be fitted by a ln(t), the short time- 
scale always shows the square-root law). Such a crossover between square root and 
exponential relaxation, occurs in all the experiments of Chiorescu et al. [85] in both 
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longitudinal and transverese fields. Note that it was also observed in the low field and 
high temperature regime (see Thomas and Barbara in [57]). 

2.2. Basic model. 

To proceed with some theoretical interpretations and conclusions, first of all, one 
should establish a basic model describing the Mni 2 spin structure. Early experimental 
works of Sessoli et al. [2,13] have reported an observation of the S=10 collective spin 
ground state, and high frequency-EPR experiments of Barra et al. [5] have allowed to fit 
its results to a "giant spin" S=10 Hamiltonian which includes fourth-order anisotropy in 
the following form: 

H G = -DS Z 2 - K,|S Z 4 + K ± (S + 4 + S. 4 ) - gu^HS, (1) 
neglecting higher-order terms, with D/k B ~0.56 K, K|/k B ~ l.llxlO" 3 K and 
K ± /k B = 2.9xlO" 5 K. Recently, a new experimental technique, the "submillimeter 
spectroscopy", was applied to the Mni 2 magnetic clusters (see [14]). This technique can 
measure directly the energy levels which, in principle, give full information about the 
Hamiltonian of the system. The information can be used to tune the Hamiltonian Eq. (1) 
with respect to higher-order anisotropy terms. 

In this model, described by Eq. (1), the fourth-order terms contribute to tunneling 
and, therefore, play a crucial role (discussion on the importance of neglected in Eq. (1) 
higher-order terms see in Section 3.2.). However, before using this Hamiltonian, it is 
necessary to understand how the model of collective spin S=10 ground state is stable 
with the increase of the temperature. In order to succeed in this task, one can try to 
calculate the energy levels of the entire problem (with all the couplings between 12 Mn 
ions as it is shown in Figurel). The dimension of the Hilbert space of such a problem is 
10 and, due to the obvious uncertainty in experimental determination of the anisotropy 
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constants, any such attempts will be rather useless. To achieve a solution of the problem 
at least in the region of the temperatures below 150 K, Tupitsyn et al. [15] suggested to 
apply an idea of the "reduced" Hamiltonian. It was assumed that the largest coupling 
constant Ji (see Figure 1), which is of the order of 200 K, locks two Mn ions with Si=2 
and S2=3/2 into a spin state of S i 2 = 1/2 up to the temperatures of the order of Ji (the 
similar tactics - "dimerization" was exploited in [16,17], followed by suggestion in 
[13]). Using this assumption, one can apply an effective 8-spin Hamiltonian Eq. (2) to 
check the temperature stability of the S=10 ground state. 

H= ^L (C 1 S H Ov+C||S^o v z )+C 2 £ a„a p ] + gn B H(£ S^ + £ o v ), (2) 

^ <jx,v> <a,p> ix=\ v=l 

where S and a distinguish four S=2 spins from the outer shell and four combined o=l/2 
spins (Figure 8). The Hilbert space of this Hamiltonian is 10 4 and now it is possible to 
use the method of the exact diagonalization to get the structure of the energy levels. In 
spite of the fact that this Hamiltonian is "reduced" (or "truncated") it includes both 
exchange and anisotropy due to the second term with Cy and simulates the whole Mni 2 
molecule up to the temperatures about 150 K at least. Note that Eq. (2) ignores 
Dzyaloshinskii-Moria (DM) interactions which could also be important in Mni 2 (see, for 
example, [18]) due to the fact that, in antiferromagnetic systems with strong couplings, 
this term can produce a single-ion anisotropy which can affect the tunneling (but this 
type of anisotropy can be "simulated" by the exchange anisotropy). In general, this is a 
time-reversal symmetry breaking term which can remove the Kramer's degeneracy of 
the energy levels (if any). Recent theoretical works (see [16] and references therein) 
suggest rather strong DM interaction in the Mn i2 clusters. 

Since we do not know the constants d, C 2 and Q, we should use some 
experimental results to determine them. First of all, we have calculated the 
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magnetization M Z (T,H Z ) and the susceptibility %n(T,H z ) for different T and H z . 
Comparing these calculated values with the experimental ones, we have tuned our 
coupling constants to get an agreement with the experiment. As can be seen from Figure 
9, we found very good agreement for Ci=-85 K, C2=55 K and Qp-7.5 K. After that we 
have used these constants to predict the magnetization and the susceptibility in a 
transverse field, M X (T,H X ) and %±(T,H X ). Figure 10 shows a good agreement with the 
experiment up to 150 K. Above this temperature our model is no longer valid since the 
pairs of spins S=2 and S=3/2 become unlocked. These calculations allow in particular to 
conclude about the temperature range of validity of the collective spin S=10 model. As 
can be seen from Figure 11, the multiplet S=9 becomes occupied above approximately 
40 K and the "giant spin" model with S=10 becomes invalid. Another calculation made 
in the simple limit of a spin S = 10 allowed to fit very well the magnetization curves, 
but only up to 30-40 K, showing also that the giant spin S=10 cannot be valid above this 
temperature [18, 21] ; interestingly, the measured (and fitted) magnetization curves 
could not be distinguished from an hyperbolic tangent with S=10, showing the Ising- 
like character of the ground state S=10, which is occupied at temperatures below 10 K. 
This does not mean that the upper levels are not occupied (the phenomenon of thermally 
activated tunneling developed below, will show that), but simply that the weight of the 
ground state S=10, is dominant on the magnetization curves below 10 K. We can 
conclude that the tunneling in the S=10 multiplet at the temperatures above 30-^40 K 
should be faster than in the S=9 multiplet (for example) since higher multiplets have 
broader and higher energy barriers. Experimentally this conclusion is confirmed by the 
observation of well-defined and equally spaced resonances at the temperatures above 30 
K [18, 21] (instead of randomly spaced resonances). It is also interesting to note that 
neutron scattering experiments performed on Mni2-ac reported the same temperature of 
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about 40 K [19] for the transition from the S=10 ground state to S=9. This means that 
the model of "reduced" Hamiltonian works well in the region of temperatures not higher 
than 120-^150 K. In the frameworks of almost the same model it was also shown by 
Zvezdin et al. in [20] that the susceptibility measured along a transverse field %j_(T,H x ), 
exhibits a peak in a transverse magnetic field at about 7t8 T. This peak can be 
interpreted in terms of a "resonance" between the two states symmetrical with respect of 
the applied field. Experimentally this transition can easily be hidden by the effect of 
fourth order anisotropy terms on the magnetization curve, and in order to avoid thermal 
resonance, the temperature should be at least 0.2 K. 

Now, after the temperature range of validity of the "giant spin" model is 
established, it is easy to estimate from Eq. (1) the critical values of the longitudinal 
magnetic field H n at which the intersection of energy levels occurs. The condition forthe 
intersection of the two levels S z =m>0 and S z =(n-m)<0 [9,10]), simply reads as ([4,11]): 

H n ~nD/g(i B . (3) 
(For simplicity's sake, we neglect here all the other terms in Eq. (1) but one can take 
them into account in a more detailled analytical expression or numerically.) Note that, 
according to the fact that the steps in the hysteresis loop were discovered only when the 
field is increased in the direction opposite to the magnetization, we have changed a sign 
before the H z term. When the field is decreased, being parallel to the magnetization, 
there is no steps because the system is in its true ground state and there is no possibility 
to tunnel to another well until the field passes through zero. The value D/k B ~0.56 K 
gives with Eq. (3) H n =0.42n T, whereas the experimental value is H n ~ 0.44n T, as can 
be seen from Figure 5. At these values of the magnetic field H z the levels m>0 and (n- 
m)<0 come into resonance and the tunneling channels open. 
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The experimental measurements of the "effective barrier" height [21] allow us to 
identify the lowest energy level (S z =m t ) where tunneling is fast enough to be recorded in 
a magnetization experiment. All the levels above this level have larger tunnel splitting 
and therefore tunnel much more rapidly. They effectively short-circuit the top of the 
barrier. This means that the "effective" height of the barrier in zero magnetic field is 
E eff (0)=D(S 2 -m t 2 ) while in non-zero longitudinal magnetic field we have E eff (H)=- 

E e ff(0)(l+H z /2DS) 2 (we neglected all other terms in Eq. (1) and assumed S»l). In the 
high temperature regime (T~2.6-^3 K), the relaxation time behaves roughly 
exponentially and follows the Arrhenius law, with E eff (H)=Tln(x(T,H)/xo). Measuring 
the relative size of dips on the curve E eff (H), one can estimate m t . As can be seen from 
Figure 12 [21], the height of the barrier in dips is reduced by about 10%, corresponding 
to tunneling from the levels m t = 3 -f 4. 

However if H z =0, all the levels are in resonance, and tunneling can in principle 
take place simultaneously from all the thermally excited levels. The tunneling rate 
r (m) between levels m and -m at zero temperature (and with no bias field between the 
levels) is approximately given by 

r (m)~A m 2 /G (4) 
(see, also [22]) where A m is a tunneling splitting, Go is a level-broadening (we assume 
Go>A m ). At nonzero temperature one has to include the thermal population of the 
excited levels, which gives for the tunneling rate 

T _1 ~£ r (m)exp(-E m /k B T), (5) 

m 

where E m is the energy of level with S z =m. Without a transverse field in Eq. (5) there 
are only five terms which correspond to the tunneling between the levels linked by 
fourth-order term in Eq. (1) (i.e., only even values of m). From this "toy" model we can 
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estimate the crossover temperature T c at which the thermal activation over the barrier is 
replaced by the tunneling through the bottom of the barrier. At T < T c the function F m =- 
E m /k B T+21n(A m ) has a maximum at m=10 while at T>T C the maximum is at m=2. Then, 
for the temperature of the crossover between these two regimes we can take T C ~(E 2 - 
Eio)/21n(A 2 /Aio). Using the method of the exact diagonalization, we have calculated 
from Eq. (1) the energy spectrum and the values of tunneling splitting A m> - m . Assuming 
that Go is independent of m (actually there is a weak dependence but we neglected it in 
this "toy" model), we get the crossover temperature T c ~1.3 K for H z =0. 

This simplified model is far from reality since it ignores correct description of the 
interactions with environment and involves only transitions Am=±4 while 
experimentally almost all the transitions with Am=± 1 are observed (see Thomas et al. 
in [11]). Therefore, this model cannot be used to explain e.g. the field dependence of the 
relaxation time. However, even from this model, it is clear that the tunneling between 
the lowest levels at the temperatures about 2 K is already unfavorable. Therefore, at this 
temperature, the relaxation process should involve at least three steps: 1) thermal 
activation (by phonons) to excited levels (for example S— >m t ); 2) tunneling across the 
barrier (m t — >-m t ); 3) transition to the true ground state with phonon emission (-m t — >- 
S). This is Thermally Assisted QTM suggested for the first time in [9,10]. 

This was in zero transverse field. In the presence of a transverse field the situation 
is completely different. The splitting A ra= s is increasing proportionally to the power 
2S=20 of the ratio transverse to anisotropy field (see [71(a)], for example), the 
relaxation should speed up very rapidly as soon as the transverse field is a sizable 
fraction of the anisotropy field. Results of Barbara et al. [18] strongly suggested ground 
state tunneling in Mni 2 if the transverse field reaches 3-4 T. Chiorescu et al. [85] 
showed that in such a field, the relaxation is fast enough to be easily and completely 
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measured. Furthermore this relaxation results from tunneling through the barrier 
between the ground states S z =± 10 and becomes faster than the relaxation by thermal 
activation above the barrier at the crossover temperature T c ~0.8 K. In this case where 
the transverse component of the applied field, is much larger than all the other 
transverse matrix elements, Eq. (1), in principle, can give quite satisfactory explanations 
of magnetic relaxation. 

This is no longer the case in low transverse fields because then, even small 
transverse matrix elements may be relevant, in particular those resulting from the 
environment which is not taken into account in this section (see section 4). 
Nevertheless, using this equation one can still predict very interesting effects which can 
be experimentally observed at low temperatures and low transverse field. The tunneling 
splitting depends on the Haldane topological phase [32] originating from the quantum 
interference of possible paths (around the hard axis) between two potential minima (S 
and -S). This topological phase can be changed by an external magnetic field, causing 
the oscillations of the tunneling splitting [33,36,24,23,31]. These oscillations were 
already experimentally observed in the similar system Feg [25] but we will discuss this 
system later. The easiest way to see the oscillations of the tunneling splitting 
analytically is to truncate the Hamiltonian in Eq. (1) to a simple low-energy 2 - level 
Hamiltonian. However, due to the presence of fourth-order terms, this task becomes 
rather complicated and to our knowledge the form of such a truncated Hamiltonian is 
not established yet. Moreover, due to the fact that higher-order anisotropy terms (up to 
20th order with S=10) can contribute importantly to the value of the tunneling splitting 
(see discussion in section 3.2), any attempt to calculate the tunneling splitting precisely 
become rather pointless because higher-order anisotropy terms cannot be measured with 
present experimental techniques). Nevertheless, since we just want to show the principle 
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things, we proceed with the simple biaxial Hamiltonian which includes an easy 
axis/easy plane anisotropy in the following way: 

H=-DS z 2 +ES x 2 -g^ B H x S x (6) 
The calculation of tunneling splitting for an isolated tunneling spin in the instanton 
technique began with the work of Enz and Shilling [71b]. For a spin tunneling in 
interaction with a background spins, results were first obtained by Prokof ev and Stamp 
[28], with much more detailed work appearing later (see Tupitsyn et al. in [23,31]). For 
different aspects of this problem see also [27,34-36,71,28,72]. The obtained 2-level 
effective Hamiltonian reads simply as: 

H eff =2A T x cos[7uS-¥], (7) 
where A is the tunneling splitting in zero external field (see, for example, 
[29,30,23,31]), x x is the Pauli matrix, and *P is the Haldane topological phase: 

^=7lg^l B H x /2[E(E+D)] 1/2 . (8) 
This expression of *P was given by Garg in [24] for the particular Hamiltonian where 
the quantization z-axis is choosen along the hard axis and the field is applied along this 
axis. Note that the result relative to the Haldane phase in [23] was given in the limit 
IEI»IDI i.e., ^TtgileHx^E. However, this limit does not affect the physics of the 
problem in general. For integer S, we get from Eq. (7): 

A H =(ilH eff l T)=A lcos( v i')l (9) 

whereas for half-integer S 

A H =A lsin( v I')l (10) 
Eq. (9) and (10) clearly show oscillations of the tunnel splitting as a function of the 
transverse magnetic field, together with the parity effect (see [35-38,23,31]). This last 
effect tells that half-integer spin does not tunnel in zero transverse field. Note that a 
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magnetic field along the easy axis (or along the medium axis which is y-axis in this 
case) does not produce any oscillations and this can be seen mathematically from the 
fact that the effective tunneling splitting A H in this case behaves like Aolcosh( v i')l (for 
details see Tupitsyn et al. [23,31])- Eq. (7) gives the general effective 2-level 
Hamiltonian describing tunneling with a transverse magnetic field ( in thiscase *P 
= v i , H+i v i , M is complex ; the contributions from hard and medium axes are denoted by *Ph 
and respectively). 

In the case of Mni 2 , where the lowest-order transverse anisotropy term is of 
fourth-order (higher order terms have not yet been determined), we can write instead of 
Eq. (8), (for H y =0): 

V=7igjj, B H 3I /T x (D,K,i,K ij S), (11) 
where T x (D,Kh,Kj_,S), the period of the oscillations along the x-axis, can be calculated 
numerically from Eq. (1). Note that this equation has two hard axes (x and y) which are 
equivalent. This means that oscillations with the same period should be seen along both 
directions. Using the exact diagonalization method, we have calculated the tunnel 
splitting A m> _ m for different m and the results can be seen in Figure 13. First of all, it is 
easy to understand that the tunneling splitting has a non-zero value in zero transverse 
magnetic field, only for levels with even values of m (which is related to the fourth- 
order anisotropy term). For all the other levels (with odd m) the tunneling splitting is 
non-zero only if the magnetic field is finite (which also produces transitions due to S+ 
and S. ). In order to see easily why the oscillations can be seen only in a finite region of 
transverse magnetic field (from -H c to H c ), let us again forget for a moment about the 
fourth order term and return to Eq. (6). Combining two nondiagonal terms, we get the 
function A(0,(p)=(sin(0)cos((p)-H x /2ES) 2 . When H X <H C =2ES, A(0,(p) as a function of 
(p has two local minima at nonzero cp. Since the Haldane phase is nothing else but the 
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area on unit sphere enclosed by two possible paths between two minima S Z =S and S z =- 
S, in the case of H X <H C the topological phase (which is an imaginary part of the 
instanton action) has nonzero value and the changes in H x results in oscillations of 
A. However for H X >H C =2ES, the function A(0,cp) has only one local minimum at cp=0, 
meaning that both paths joining the states S and -S coincide (up to the cp-fluctuations of 
trajectories which renormalize the value of H c to [2E(E+D)] 1/2 S). The area enclosed by 
these two paths, i.e. the imaginary part of the instanton action is zero and, therefore, 
there are no more oscillations of the tunneling splitting (see also [24]). 

It is important to note that the number of zeros of the A m> . m (Hj_) function strongly 
depends on the symmetry of the anisotropy terms. For the Hamiltonian of Eq. (1) (case 
of Mni 2 ) where the lowest transverse anisotropy term is of fourth order, the number of 
zeros for even values of m (along positive or negative direction of field) should not 
exceed the number of times (v=l,...,5) the operator S. 4 has to be applied to the state 
S z =m to reach the state S z =-m. (However, this is only valid for K_|_>0 ; if Kj_<0, the axes 
x and y are no longer the hard axes and A(Hj_) obviously would not show the 
oscillations along these directions.) The tunneling splitting for odd m has an additional 
zero at zero magnetic field. In this case the chain of operators S_ 4 or S+ 4 which should be 
applied on the state Im) to reach the state l-m) must be completed by additional 

operators S- (or S+) which come from magnetic field term (i.e. l-m) =(S.) 2 (S. 4 ) v lm> with 

V=0,...,4). The same situation should occur for the tunneling between the levels 
involved into the resonance by applying non zero longitudinal magnetic field (say, 
levels S z =-m and S z =m-n). Only the levels linked by (S + 4 ) v or (S. 4 ) v (v=0,...,4) can have 
non zero tunneling splitting in zero transverse field (see Figure 14). Moreover, it is easy 
to see from this figure how the period of the first oscillation decreases when n increases 
inside the group of curves (-10, 9) -K-10, 7) (or (-10, 5) -K-10, 3)). To reach the states 
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S z =9, 8, 7 from the state S z =-10 we should apply S. 4 four times and then, to complete 
the chain linking the mentioned states, apply S. three, two or one time. 

In addition, if we apply the transverse magnetic field at different azimuth angles (p 
and increase (p from 0, the amplitude of the oscillations becomes smaller and vanishes at 
cp=7l/4 since in the case of fourth-order anisotropy 7l/4-axis is an easy axis. However, 
increasing cp from 7l/4 to Till, one can see the same curves. This is quite obvious from 
the symmetry of the problem (see Figure 15). This is the case when Kj_>0. If K_|_<0, as 
we noted already, x and y axes are no longer the hard ones. Two axes along the 
directions of cp= ± 7t/4 become harder for the system. Some oscillations along these 
directions should be observed. 
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3. Feg octanuclear iron (III) complexes. 
3.1. Experimental results. 

Another molecule, so-called Feg, with the chemical formula 
[Fe 8 02(OH)i 2 (tacn) 6 ] 8+ , where tacn represents the organic ligand triazacyclononane, is 
under intensive investigations at the present time. It contains eight iron (III) ions 
(S=5/2) with strong antiferromagnetic couplings between ions [39] (see Figure 16). 
Similarly to Mni 2 they form an uncompensated S=10 collective ground state. Four Feg 
ions in the middle of the molecule are in the so-called "butterfly arrangement". This 
system is nearly orthorombic with a strong Ising-like anisotropy and giving an energy 
barrier of about 24 K (i.e., about one-third of that in Mni 2 ) [39,40]. An analysis of the 
magnetic susceptibility v.s. temperature shows that only the levels S>8 are populated 
near 10 K [39]. (As for Mni 2 (see above), one must say that upper levels are also 
occupied, but the weight of the ground state S=10, dominates the susceptibility.) 

Magnetic relaxation experiments have been performed in Feg, following the same 
procedure shown above, for Mni 2 . The relaxation rate becomes temperature- 
independent below 0.35 K, as was shown by Sangregorio et al. in [41]. This can also be 
seen in the hysteresis loop recently shown for a single crystal [42] (Figure 17). Like in 
Mni 2 , equally-spaced steps were observed, but with a smaller spacing : AH -0.22 T 
instead of 0.44 T. In [41] it was found that at these values of field (H n = 0.22n T), the 
relaxation becomes much faster than that in the plateaus. These observations give a 
second example of tunneling across the anisotropy energy barrier, when the levels from 
the opposite sides of the barrier come into resonance. At low temperature and without 
transverse field, the ratio of the relaxation time measured at resonance and out of 
resonance is larger than in Mni 2 by one or two decades, showing that the relaxation in 
Fe 8 is, in this case, faster than in Mn i2 by the same factor (one or two decades). 
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(However, as mentioned above, in the presence of a transverse field of a few Tesla the 
relaxation of Mni 2 becomes much faster ; if the transverse field is of 3-4 T, ground state 
tunneling between S=+10 and S=-10 is observed below the crossover temperature 
T c ~0.8 K [85].) In Fe 8 , ground state tunneling occurs below T c ~0.35 K, but in this case 
it is not necessary to apply a transverse field. To fit the relaxation data, a stretched 
exponential law M(t)=M(0)exp[-(t/x) p(T) ] (with (3(T) increasing from approximately 
0.4-h0.5 (below 0.4 K) to nearly 1 at T~l K) was used [41]. This law was also observed 
for measurements on an oriented crystal by Ohm et al. (see [43,56]). At temperatures 
below 0.35 K and at short times, the best fit of the data is a square root behavior 
M(t)=M(0)[l-(t/x S hort)] 1/2 ) (see Figure 18). This law was, in fact, predicted by the theory 
of Prokof'ev and Stamp [44] for the relaxation due to tunneling at the bottom of the 
barrier, at short times and low temperatures, with initial magnetization near saturation. 
Later these measurements were repeated at lower temperatures (T=40 mK) by 
Wernsdorfer et al. and the square root law was confirmed [45]. Note, however, that this 
law was also found in a zero-field cooled annealed sample, which is then allowed to 
relax in a finite field, i.e. near zero magnetization. A similar effect is also observed in 
Mni2 [85] (see below). This needs further theoretical investigations (see [91]). 

In conclusion, in Fe 8 the relaxation goes in the pure quantum regime via the 
tunneling through the barrier between the states S z =± 10 at temperatures T<T C ~0.35 K 
and zero transverse field. Similarly, in Mni 2 the ground state tunneling was also 
observed, however a transverse field of 3-4 T must be applied to achieve this result ; the 
crossover temperature is then larger than in Fe 8; (T c ~0.8 K for Mni 2 ). The application of 
such a transverse field in Fe 8 would makes the relaxation so fast that it would be 
impossible to measure it (unless by EPR). At higher temperatures ac-susceptibility 
experiments of Caneschi et al. [42] show peaks similar to those observed in Mni 2 
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[12,21], but they are more pronounced than in Mn i2 , from which we may conclude that 
m t (m-value for the barrier short-cut) is larger in Feg. However one should keep in mind 
that the collective spin S=10 of this molecule breaks down near 10 K, suggesting that 
the peaks observed at high temperature (7 K) come, for a large part, from multi-spins 
tunneling. This is corroborated by the fact that the peaks are not regularly separated and 
that their mean separation looks closer to 0.14 T than to the 0.22 T of the spin S=10. 

3.2. Basic model. 

High-Frequency EPR technique was applied by Barra et al. [40] to investigate the 
magnetic anisotropy in Feg molecules. They found a biaxial anisotropy described by the 
Hamiltonian H=-DS z 2 +E(S x 2 -S y 2 ). More recent neutron spectroscopy experiments (see 
Caciuffo et al. in [46]) have reported the presence of fourth-order term. According to 
these experimental data, the Hamiltonian for the Fe 8 molecule can be written as follows: 
H G = -D S Z 2 + E (S x 2 -S y 2 ) + K ± (S + 4 + S. 4 ) - gu*HS (11) 
which is equivalent to: 

H G = -DS Z 2 + ES X 2 +K ± (S + 4 + S. 4 ) - gu*HS (12) 

where D/k B =(D -E )/k B ~0.23 K , E/k B =2Eo/k B ~ 0.092 K and K ± /k B --2.9xlO" 5 K (g ~ 2, 
see [40], [46]). Eq. (11) and (12) are valid only in the temperature range where the 
collective spin S=10 can be defined, i.e. at temperatures lower than 10 K [39] (this 
temperature is 40 K in Mni 2 ). 

Similarly to the case of the Mni 2 molecule we can estimate the value of the 
longitudinal magnetic fields at which the levels from the opposite sides of the barrier 
(say, S z =m and S z =n-m) come into resonance. Eq. (12) gives: 

H n =(nD/g^i B )[l+E/2D]. (13) 
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(we again neglected the fourth-order term). The contribution E/2D compensates the 
difference between D and D (compare with Eq. (3)) and therefore one must have 
D =D[1+E/2D], which is true and identical to D=D -E (see above). Eq. (13) gives 
H n ~0.205n (in Tesla), whereas the experimental value is 0.22n T. The above 
considerations show that we may have some uncertainty in both the actual values of the 
constants and in the types of the anisotropy terms included into the Hamiltonians. As 
was pointed out in [44,47], higher-order transverse anisotropy terms (even with very 
small constants) can make an important contribution to the value of the tunneling 
splitting. This can be seen easily from the perturbation theory (for the lowest order 
perturbation approach for the tunneling splitting see [48,49,71]). Following [71] a 
simple form for the tunnel splitting can be written [21] (we omit here dependence on 
value of S): A„ m , m ~ D(K_|_ P /D) 2m/P (where p is the order of the anisotropy term in the 
Hamiltonian). As an example, Eq. (12) gives for m=S=10 and p=2 or p=4, A_i ,io ~ 
D(E/D) 10 or Ai ,-io ~ D(K ± /D) 5 . The contribution of e.g. the 10th order term gives 
already D(Kj_io/D) . Unless the quasi-exponential increase of A with p, a divergence of 
higher orders of the tunneling splitting is forbidden due to the fast decrease of the 
constants Kj_ p (these constants are directly connected to the crystal field parameters 
which we know, decrease very rapidly with the expansion order p). Nevertheless, A 
depends on the value of Kj_ p in a so crucial way that all the terms up to 20th order can be 
important. However, the parameters Kj_ p of the higher orders (except p=4) are 
experimentally unmeasurable at the present time. This makes very problematic any 
quantitative calculations of the tunneling splitting from Eq. (12). The actual behavior of 
A m , n -m(H x ) strongly depends on the values of the anisotropy constant. Any uncertainties 
in these constants cause the changes in the period and amplitude of the oscillations in 
the transverse magnetic field (see Figures 13-15). Furthermore, number of the 
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oscillations (which are confined in a given field interval) itself depends on the relative 
values of the constants. In the case of the simplest example of the biaxial anisotropy of 
Eq. (6) the field interval for oscillations is [-H C ,H C ], where H C ~[2E(E+D)] 1/2 . Consider 
the tunnel state where 1-10) and 110) are admixed, with splitting Ai ,-io- In order to 
reach the state 1-10) starting from the state 110) , the operator S+ p must be applied 20/p 
times (for even p). The main consequence, is that the number of oscillations depend 
directly on the symmetry of the anisotropy: v(p=2)=10 oscillations for a second order 
anisotropy or v(p=4)=5 oscillations for a fourth order one. It is clear that if Ki increases 
from zero, the transition from v(p=2)=10 to v(p=4)=5 will not be discontinuous, 
because both periodicity will be involved in the interference. If Kj_ is negligible (in Eq. 
12), only the second order term contributes to the splitting. However, as Kj_ increases 
above some critical value K_|_ C ~2E/S 2 (with Kj_>0) its contribution becomes dominant 
and v(p=2)=10 decreases progressively to the value v(p=2)=4 (as in the case of Eq. (1)). 
Simultaneously, the period and the amplitude of oscillations increase. One must note 
that in all that the changes in anisotropy constants are not small: the transition between 
v(p=2) and v(p=4) occurs when the parameters of consecutive orders in Eq. (12), are 
nearly the same. In the above example, K_|_ C ~2E/S 2 -1.8 10" 3 K is about 10 2 times the 
real value of Kj_! Unless very unusual values of the crystal field parameters, such a 
transition could not be observed. In any case, this situation is not stable since the 
splitting resulting from such large Ki is itself very large: A c ~ 7 10~ 3 K, i.e. nine orders 
of magnitude larger than the actual splitting (given above; let us note casually the huge 
effect of Kj_ on A). Comparing the value A c to the splitting given e.g. Figure 22, it is not 
difficult to imagine that the energy spectrum will changes dramatically and all the levels 
with definite S z will be completely admixed (S z will not be conserved). In the opposite 
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case (Kj_<0), the increase of IKJ above Kj_ c leads to a disappearing of the oscillations 
because the x-axis is no longer the hard axis of the system (as well for the y-axis, due to 
the tetragonal symmetry of Eq. (12)). 

Any other anisotropy terms of order higher than four, could in principle, change 
the number, the period and the amplitude of the oscillations. These oscillations have 
been recently observed in the system Fe 8 by Wernsdorfer et al. [25]. Figure 19 shows 
the measured tunneling splitting as a function of transverse magnetic field at different 
azimuth angles cp. Using the Hamiltonian (Eq. (12)), we have calculated the tunnel 
splitting in a transverse field by simple diagonalisation of the 21x21 matrix. We get 
similar oscillations as in [25] with the following parameters: D/k B =0.23 K , E/k B =0.094 
K and K ± /k B =-3.28xlO" 5 K. There are some differences between the curves calculated 
for cp=0° and the measured one (the calculated curve is sharper near the nodes, and the 
experimental curve shows some increase of the value of A_io,io(H x ) in the nodes). 
However, the curve calculated for cp=l° (shown in Figure 20), is more similar to the 
experimental one. This suggests some misorientations in the experiments. The reason 
for this suggestion is connected with the fact that an increase of the magnetic field along 
the medium axis (which is the y-axis in our case) produces an increase of the tunneling 
splitting which obeys the relation Aolcosh($ y )l (where 3> y is the Haldane phase for the 
magnetic field H y and A is the tunneling splitting in zero field). The simplest way to 
simulate this non-zero component of H y field (not unique however) is to introduce some 
"misalignment" angle between the x-axis and the direction of the applied magnetic 
field. One cause could be the mosaic which is always present in molecular crystals of 
this type, which is also of the order of 1°. However this gives random misorientations 
between say +0.5° and -0.5°. Such a mosaic should modify the curve A_i ,io(H x ), but in 
a distributed way, and this could contribute to the observed broadening of the nodes. 

23 



Another origin for misorientation comes from the triclinic symmetry of the Feg 
molecule; in this case symmetry lowering should induce new crystal field parameters 
and new contributions to the tunneling splitting, which are in first approximation taken 
into account by field H y . In what follows, we will use the value of this "misalignment" 
angle m =l°. 

The calculated A_i ,io(H x ) with zero longitudinal field ( but m =l°), can be seen in 
Figure 20. The period of the oscillations is about 0.41 T in agreement with the 
experimental value. The curves with the larger value of cp (up to %I2) clearly show that 
the oscillations disappear (in agreement with the experimental behavior) when the 
direction of the applied field approaches the medium axis (y). As for the absolute value 
of A_io,io(H x ) and the shape of the last oscillation, the agreement is not very good. These 
discrepancies may give an indication that the Hamiltonian Eq. (12) is not quite 
satisfactory with respect to the unknown higher anisotropy terms. It would always be 
possible to choose the value of these terms up to highest order (20) to get the best fit. 
One might also consider the effects of couplings to the environment. In order to give an 
idea of the influence of high order terms, we show Figure 21, the effect of the fourth- 
order term added to the second order one. Starting from Eq. (12), we have calculated 
numerically the period of oscillations T H , for different values of Kj_. When this term is 
null i.e. when the period is given by Eq. (8), the value of T H is about one half the 
measured one. The value of the fourth order term Kj_ =-2.9 10" K, allows to recover the 
measured period. One could also take another value of Kj_, determined independently, 
and fit the period on the eight order term (for example). This is a kind of unproductive 
activity (at the present time), and we stop here this discussion. 

At higher values of Ki the period depends on Kj_ almost linearly but the value of 
interest (IKJ~3-h4xlO" 5 K ) is in the nonlinear region. In the absence of longitudinal 
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field, (A_ m>m . n with m=10 and n=0), it is easy to interpolate T H by some simple formula 
even in this region, but all the other anisotropy constants also have uncertainty in their 
values. To tabulate T H as a function of all available anisotropy constants, we have used 
the combination of the method of exact diagonalization with the method of polynomial 
interpolations. The final formulas are valid in quite a wide region of the values of the 
anisotropy constants (in Kelvins): 

D/k B e [-0.06,-0.45]; E/k B e [0.05,0.13]; K ± /k B e [-0.8xl0" 5 ,-5.2xl0" 5 ] (14) 
and can be written as follows: 

T H =(2k B /g^i B )[E(E+D)] 1/2 £ X, K G^ V X 

IX,V=l 
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For the general case (n>0), one has to come back to full numerical calculations. 
The period depends on the values of n and m (Figure 22 and 23 a,b). Figure 22 shows 
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the tunneling splitting A m _ m (in Kelvins) as function of H x (in Tesla) in zero 
longitudinal magnetic field for the azimuth angle (p=0 (with "misalignment" angle 
m =l ). Due to the presence of the second-order term in Eq. (12) all the A m _ m have 
nonzero values in zero transverse field. As we discussed above, negative Ki tends to 
decrease the number of the oscillations along the hard x-axis because it makes x- 
direction easier for quasiclassical motion of the giant spin. This effect leads to the 
decrease of the imaginary part of the instanton action (by decreasing the area on the unit 
sphere, enclosed by two possible paths joining quasiclassical minima - S and S). In this 
case it is difficult to calculate the exact number of oscillation. First of all, an analytical 
solution of the problem in the presence of the fourth-order term is not easy, in particular 
the perturbation theory cannot be applied in large transverse magnetic field. The 
instanton technique and WKB approximations give non-analytical solutions which 
require numerical calculations (there is no reason to apply it since we already have an 
answer given by the method of exact diagonalization). As mentioned above, the only 
conclusion that we can make here is that the number of oscillations is determined by the 
exponents oc,p\y in the chain of operators (S+ 1 ) a (S± 2 ) |3 (S± 4 ) T which should be applied to 
the state lm> to reach the state In-m). The operator S+ comes from the transverse 

magnetic field term, S+ 2 and S+ 4 come from the second and fourth order anisotropy 
terms, respectively. It is clear that the whole picture is defined by the combined 
symmetries of the anisotropy terms. This is apparent in Figure 22, with the change of 
the period with m and the shape of the last oscillation. Note in particular that the 
quantities Ai 0> -io, A 8j _ 8 and A 6; _ 6 (with the difference 8S Z =4 in the lengths of the chains 
connecting lm> and l-m) states) show a similar behavior. The same situation takes 
place for Ag r g, A 7j . 7 and As,^. For smaller m, the structure of interference is also affected 
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by the admixing of states with different m (S z is not a good quantum number for the 
Hamiltonian Eq. (12)). 

The same conclusions can be drawn, in the presence of a longitudinal field. The 

functions A_io,io-n(H x ) calculated for cp=0, are plotted in Figure 23 (with n even - (a) or 

odd - (b)). The resonant longitudinal field for particular n is defined by Eq. (13) (of 
course, it was necessary to tune the value of the field around these H n ). The levels S z =- 
10 and S z =10-n with odd values of n cannot be linked in zero transverse field (there are 
no matrix elements between them) and therefore they have a zero tunneling splitting at 
H x =0. The levels S z =-10 and S z =10-n with even values of n can be linked even in zero 
transverse field, and they have a finite tunnel splitting. (As an example for n=l and n=2, 

one has to apply S+ and S+ 2 respectively.) The curves calculated in Figure 23, can be 

compared with the ones measured in Fe 8 [25] and plotted in Figure 24. The absolute 
value of the tunneling splitting, again, is not in the best possible agreement with the 
experiment. Unfortunately, the authors of [25] did not show behavior of the curves with 
n=l,2 (and with larger n) at higher values of transverse magnetic fields. The comparison 
between measured and calculated curves in the region of high fields can help to 
determine which types of anisotropy are important in the "giant spin" Hamiltonian. In 
other words, measurements of the tunnel splitting as a function of a magnetic field can 
be used to determine the crystal field parameters of the "giant spin" Hamiltonian. This 
method should be quite sensitive to the anisotropy of the higher orders 

In addition, we would like to stress that the actual behavior of the tunneling 
splitting in the region of the nodes is very sensitive to the environment. This point will 
be discussed below in the part of this review related to the low temperature limit of the 
relaxation of the magnetization (Environmental effects). 
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Figure 1. Interaction scheme of Mni 2 molecule. 

Figure 2. Variation of x(H z ) measured in [7,9] at T>T C . Note that, in this figure (which 

we take from [9]), the value of Tln(x(H z )/x ) is plotted to show also that Uo(H z ) 

(effective barrier height) deviates from ~DS 2 (1+H Z /2DS) 2 . 

Figure 3. Hysteresis loops of Mni 2 with the field along z-axis from [11]. 

Figure 4. Field variation of the derivative 3M Z /3H Z taken at 1.9 K along the hysteresis 

loop of the single monocrystal [12]. The sharp peaks correspond to the magnetization 

jumps and flat regions correspond to the plateau located between the jumps. The 

continuous line is a fit to the Lorentzian peaks centered at H n . 

Figure 5. Relaxation times at temperature 1.9 K versus H z obtained in [11] from 
repeated measurements for given H z and T on the hysteresis loop. The insert shows the 
relaxation time drops against inverse temperature. 

Figure 6. Field dependence of the relaxation time x=%"(co) /co(x'(co) - %'(°°)) from AC 
susceptibility measurements [12]. The dashed line represents the fit of the thermal 
activation background and the continuous curve is a fit of thermally activated resonance 
dips. 

Figure 7. (a): Hysteresis loops obtained on a single crystal of Mni 2 (main phase) from 
torque experiments [85] performed at temperatures between 1.3 and 0.4 K. The 
magnetization was first saturated in a large positive field. The field was then decreased 
to zero and reversed. Data points were taken between -2.5 and -5.5 T. The sweeping 
field velocity was equal to 10.8 mT/s. The steps amplitude depend on temperature, but 
only above 0.5 - 0.8 K. Below this temperature the hysteresis loops are independent of 
temperature, suggesting tunneling from the ground-state S=10. (b): an example of 
magnetic relaxation experiments performed near the maximum of the first resonance in 
the presence of a transverse field of about 4 T, on a single crystal of Mni 2 . The data 
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were taken after saturation in a positive field and fast application of the field to the top 
of the resonance. This curve shows that the relaxation follows a square root law at short 
time-scale and an exponential law at long timescale. Insert: temperature dependence of 
the square root relaxation time. Below 0.8 K, the relaxation time is independent of 
temperature showing the existence of grounde-state tunneling between S=10 and S=-10 
in this bulk phase of Mni 2 . 

Figure 8. Simplified coupling scheme of Mni 2 molecule. 

Figure 9. (a) Temperature dependence of longitudinal susceptibility Triangles are 
experimental results from [59]. Dashed line is %|T calculated from Eq. (2) with the 
contribution only from the energy levels up to 100 K starting from the ground state (85 
levels). Solid line is %|T calculated from Eq. (2) with the contribution from all the 
energy levels (10 4 levels), (b) Dots are magnetization curves <M Z >/M S versus 
longitudinal field measured in [59] at different temperatures. Solid curves are the same 
curves calculated from Eq. (2). 

Figure 10. (a) Temperature dependence of transverse susceptibility Triangles are 
experimental results from [59]. Dotted line is %j_T calculated from Eq. (2) with the 
contribution only from the energy levels up to 100 K. Dashed line is the same but with 
the contribution from energy levels up to 500 K (2982 levels). Solid line is %j_T 
calculated from Eq. (2) with the contribution from the energy levels up to 1000 K (8362 
levels), (b) Dots are magnetization curves <M X >/M S versus transverse field measured in 
[59] at different temperatures. Solid curves are the same curves calculated from Eq. (2). 
Figure 11. Energy spectrum, calculated from Eq. (2) up to 180 K. The stars show a 
parabolic behavior (A(S 2 -S Z 2 ) where -S < S z < S and A=0.627 K) of S=10 multiplet. 
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Figure 12. Effective energy barrier of Mrii 2 molecule versus longitudinal field from 
[21]. The solid, dashed and dotted curves represent the law E e ff(H)=Aef^-A(l+H z /2DS) 2 
for different values of the barrier height Uo=A. 

Figure 13. Tunneling splitting A„ m , m (Mni 2 ) calculated from Eq. (1) for different m 
versus transverse magnetic field. 

Figure 14. Tunneling splitting A_i ,io-n (Mni 2 ) calculated from Eq. (1) for different n 
versus transverse magnetic field. 

Figure 15. Tunneling splitting Aio,-io (Mni 2 ) calculated from Eq. (1) for the different 
values of the azimuth angle (p versus transverse magnetic field. 
Figure 16. Interaction scheme of Fe 8 molecule. 

Figure 17. Hysteresis loops recorded on a single crystal of Fe 8 molecules in [42] at 
different temperatures and at constant sweep rate 3H/3t=0.14 T/sec. 
Figure 18. Square root of time relaxation curves for Fe 8 single crystal measured in [43] 
(see also [56]) at 150 mK for M in =M s . The insert shows the distribution of T^t" 1 
extracted from the above data as a function of field. 

Figure 19. Ground state tunneling splitting Aio,-io measured for several azimuth angles 
(p versus transverse magnetic field [25]. 

Figure 20. Ground state tunneling splitting Aio,-io calculated from Eq. (12) for several 
azimuth angles (p versus transverse magnetic field. 

Figure 21. The period of the oscillations of the ground state tunneling splitting Aio,-io 
calculated from Eq. (12) as a function of the fourth-order anisotropy constant Kj_. 
Figure 22. Tunneling splitting A„ m , m (Fe 8 ) calculated from Eq. (12) for different m 
versus transverse magnetic field. 
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Figure 23. Tunneling splitting A_io,io-n (Feg) calculated from Eq. (12) for different n 

versus transverse magnetic field, (a) - even and (b) - odd values of n. 

Figure 24. Tunneling splitting A_io,io-n (Feg) measured for n=0,l,2 [25]. 

Figure 25. Quantum hole digging in the initial distribution of internal fields [53]. Note 

that these measurements were performed on a minor species of Mn^. (For details about 

the different species of Mni 2 , see [53,42] and references therein.) 

Figure 26. The digging time dependence of the hole shape in Mni 2 from [53]. 

Figure 27. Temperature dependence of the hyperfine line-width (Oh yp ) in the crystal of 

Mn 12 [79]. 

Figure 28. Tunneling distribution in Fe 8 (which, according to [44], is proportional to 
the distribution of P(^h) of the internal bias field £ H ) for annealed sample measured in 
[45]. The insert enlarges the region of the fields around the hole. 

Figure 29. Dependence of the hole shape (in initial distribution of the dipolar fields in 
Fe 8 ) on degree of annealing [45]. At initial magnetization M in < I0.5MJ the hole becomes 
independent on future annealing and gives the line- width Ohyp~1.2-r 1.6 mT. 
Figure 30. Temperature dependence of the hyperfine line- width width (Oh yp ) in the 
crystal of Fe 8 [80]. 

Figure 31. Time decay of the magnetization in Mni 2 , plotted as a function of the square 
root of the time [58]. Lines show linear fits. 

Figure 32. Square root of time relaxation curves for Fe 8 crystal measured at 40 mK in 
[45] for M in =0. 

Figure 33. The value of the tunneling splitting As, n -s versus transverse magnetic field H x 
(cp=0°) around the first node at the different values of the initial magnetization M in [83]. 
Figure 34. The value of y calculated from the Eqs. (12) and (61) in the case of strongly 
annealed sample (the Gaussian distribution of the dipolar fields), (a) - y versus 
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transverse magnetic field at different values of Gaussian half-width E D ; (b) - enlarged 
region of transverse magnetic field around the first node; (c) - the value of y in the node 
versus E D 2 . 

Figure 35. The value of y versus transverse magnetic field calculated with small 
"misalignment" angle m =l° at different values of E D . 
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4. Environmental effects. 

At the present time, current theories are not able to explain relaxation experiments 
in all the temperature ranges. Therefore, it is reasonable to consider the low 
temperatures (ground- state tunneling regime) and the higher temperatures (thermally 
assisted tunneling regime) separately. The first is better understood by the theory of 
Prokof ev and Stamp [44,47,72,73,77] while progresses in understanding of the second 
are more related to the works of Luis et al. [50], Fort et al. [51] and Leuenberger et al. 
[52]. Since the barrier height in Fe 8 is about three times smaller than in Mn^, the values 
of the tunneling splitting are larger for all states. In particular one expects a barrier cur- 
off at larger m t in the high temperature regime, and faster relaxation in the low 
temperature regime. These circumstances make Feg molecules very attractive for the 
study of the ground-state tunneling regime (i.e. pure quantum regime). The thermally 
assisted regime is usually (historically) the main subject of investigations on Mni 2 
molecules. 

4.1. Experimental picture. 

To deal with relaxation one has to take into account the environmental effects. The 
environment of a molecule is essentially constituted of ensembles of bosons and 
fermions coupled to the spin of the molecule. In real molecular systems the most 
important environmental effects come from phonons, nuclear spins, dipolar fields. The 
environment is able to absorb finite variations of energy and angular momentum. This is 
extremely important because the non-conservation of these quantities can forbid the 
tunneling. The environment is also responsible for the broadening of resonance lines. 
The shape of these lines is usually obtained from the plot of dM z /dH z v.s. H z , as this 
was first defined in Ref. [11]; in Mni 2 nearly Lorentzian resonance line-shapes were 
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found of width (40mT to 100 mT). Furthermore, as mentioned by Barbara et al. [18,11], 
tunneling resonance could not have been detected with the poor field resolution of 
conventional SQUID magnetometers, if resonance lines were not importantly broadened 
and since pure phonons broadening is too small by a factor of 100 [65], broadening of 
magnetic origin must play the major role (the line-with in the absence of environment is 
~ tunnel splitting A) . This is true at low temperature only because the effect of phonons 
is obviously dominant in the thermally assisted regime. Very similar transitions 
broadening were observed in Fe 8 by Ohm et al.; at low temperature they are of the order 
of 12-=- 15 mT [43]. 

Now we describe the origin of the resonance line-width in molecular crystals, in 
the two limits of high and low temperatures (see also [21] for qualitative description 
from experiments). In the thermally activated regime (where A at the bottom of barrier 
can reach values of the order of 0.5 K) the line-width is at least equal to the value of the 
tunneling splitting. The resonance line in this regime must be homogeneously 
broadened with Lorentzian shape due to equilibrated spin-phonons transitions, and this 
was observed in Mni 2 (see e.g. [11, 12, 60]). Note that, deviations from a Lorentzian 
line-shape can be obtained depending on the magnetic history of the sample, and this is 
because the width of distributions of dipolar and hyperfine fields are of the same order 
with the tunneling splitting. It might be surprising to see that magnetic history, can be 
important in a basic effect such as tunneling. In fact the magnetic history is always 
important when irreversible process come to play a role, and this is the case here. It is 
well known that frozen distributions of hyperfine and dipolar fields (quenched from the 
super-paramagnetic state) have Gaussian distributions and if the temperature is low 
enough to prevent fast spin reorganizations. In this case resonance lines will be 
inhomogeneously broadened, with no Lorentzian and eventually Gaussian line-shapes. 
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In the low temperature regime, where thermal fluctuations are essentially frozen, 
tunneling will only be possible through internal fields fluctuations, of amplitude H f . 
These fluctuations are due to nuclear spins in weak dipole-dipole interactions, and with 
strong hyperfine interactions with electronic spins. Internal field fluctuations involve 
spin-spin quantum dynamics of the considered system which can be tested on either 
electronic or nuclear spins. In Mni 2 each ion has a nuclear spin, whereas in Fe 8 this is 
the case for only 2% of the ions (Fe 57 ). However, each Fe 8 molecule contains 120 
hydrogen, 18 nitrogen and 8 bromine atoms, which all have nonzero nuclear magnetic 
moment. They produce (together with the dipolar moment of each molecule) the 
fluctuating internal field acting on each molecule. The amplitude H f of these fluctuating 
fields are extremely small (they can be evaluated from NMR experiments [86]; e.g. 
protons in Fe 8 give H f ~1. 4 mT [21]). Quantum tunneling being only possible within 
this range of fluctuating fields, spin reversals from S to -S will occur only in a narrow 
energy window of width H f , "digging a hole" in the initial distribution of internal fields. 
The frozen distribution of internal fields (resulting from electronic and nuclear spins 
which are not affected by the quantum dynamics) could be obtained from the measured 
dM z /dH z v.s. H z , in sweeping the longitudinal field, as we said above. However, 
contrary to the high temperatures case, the distribution which is probed here for each 
value of H z , is inhomogeneous, and it will be possible to probe this distribution as long 
as Hf will remain much smaller than its width. This will be the case unless the 
temperature increases to the point where H f becomes of the order of the total local field 
distribution. The cross-over between inhomogeneous Gaussian-like to homogeneous 
Lorentzian-like distributions is when the temperature becomes large enough to 
equilibrate the spin system. In Mn^, the tunneling window H f was evaluated assuming 
oscillations of the mean dipolar and hyperfine field with Am =AI= ± 1 (see [21,87]). 
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Recently low temperature experiments (0.04-^0.3 K) of Wernsdorfer et al. [53], 
using the "hole digging" technique [45] based on the theory of Prokof'ev and Stamp 
[44], allowed to measure the tunneling window. In this theory the relaxation rate of the 
magnetization r sqrt (H) is proportional to the distribution of internal bias field P(^h)- On 
the experiment, during the digging time t dig a small fraction of molecules (that are in 
resonance in applied external digging field H dig ) tunnel reversing the direction of their 
magnetization. This causes rapid transitions of molecules which are close to the 
resonance around H dig . Such transitions are effectively "digging a hole" in an initial 
distribution of internal fields (see Figure 25). The hole widens in time, depending on the 
digging time t d i g . The hole-width (which is obtained by the linear interpolation to t dig =0) 
gives an intrinsic broadening of the nuclear fields Oh yp (see Figure 26). It was found [53] 
that the width of the Gaussian distribution of hyperfine fields Oh yp ~12 mT which is 
larger than the value of the tunneling splitting Ai ,-io by orders of magnitude. Oh yp is the 
temperature independent up to approximately 0.4 K and then starts to increase (see 
Figure 27). (The last fact needs further theoretical investigations since the square-root 
theory [44] does not work at so high temperatures.) We should note, however, that these 
"hole digging" experiments were performed on the minority phase of Mi 2 and in the 
presence of a magnetic field this system can have a lower crossover temperature. Low 
temperature studies of main phase of Mi 2 , were recently possible. A large magnetic field 
was applied to get ride of the minority phase (see [85]). 

Similarly to M i2 it was observed in Fe 8 that the first maximum in the relaxation 
rate is not in zero applied field, but at 8 mT. As in Mi 2 this was attributed to the effect of 
internal fields [41]. This was in a powder sample. In a single crystal it was found that 
the resonance width (of about 12-^15 mT) is as in Mi 2 and for the same reason 
(intermolecular dipolar interaction), orders of magnitude larger than expected without 
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environment or with phonons only [18,11,65]. A difference between these two systems 
is that the first resonance is observed in a negative field in Mi 2 and in a positive one in 
Feg. This was interpreted as a consequence of the competition between the 
demagnetizing field -NM (shape-dependent, where N is the demagnetizing factor) and 
the local Lorentz field +(4tc/3)M [21]. Since the Mi 2 crystals are elongated, the 
demagnetizing field (which is antiparallel to M) is smaller than the Lorentz field and the 
internal field is parallel to the magnetization M: one has to apply a negative field to 
cancel the internal field. The situation is just opposite with Fe 8 . This discussion was 
relative to the most probable values of internal fields. Now, we know that internal fields 
are distributed with a width, measured in both systems (see above). Hole-digging 
experiments similar to those described above for impurity phase of Mi 2 were performed 
in Feg [45], also at low temperature. It was found that, for thermally annealed sample 
(down to values of magnetization of about -0.2M S , where M s is the saturated 
magnetization), the distribution P(^h) of the internal bias field ^ H is very accurately 
described by a Gaussian function (Figure 28). As a matter of fact such a distribution is 
expected from theoretical models, but only in the limit of high spin concentrations 
(dipolar field distribution of dense set of randomly oriented spins, see [88,54]). In the 
other limit of dilute static dipoles, the distribution has to be Lorentzian, as this was 
shown by Anderson in [55]. A study for different concentrations was done showing that 
in between these two limits more complicated and sometimes bi-modal distributions 
could occur [89]. Interestingly, in Figure 28 the maximum of P(^h) is shifted from H=0. 
The shift is even larger than in [41] and this is because the experiment was done at 
lower temperature: the magnetization and thus the local field are larger (the sample was 
not absolutely annealed) The width of this distribution (Odi p ) was found to be of the 
order of 50 mT). 
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The relaxation of the molecules, like in Mni 2 , digs a hole in P(^h) at the value of 
the applied field. The width and depth of the hole change with the waiting time i.e. with 
the time during which the field is applied (see insert in Figure 28). For initial 
magnetizations close to the saturation the hole is large and asymmetric, whereas for the 
initial magnetization less than I0.5MJ it becomes symmetric, independent on further 
annealing (on initial magnetization) and has a width of about 1.2 -r 1.6 mT (see Figure 
29). Similarly to Mni 2 this line-width is temperature independent up to 0.4 K and then it 
starts to increase (see Figure 30). Some values of the hole width were predicted by 
Prokof'ev and Stamp at zero Kelvin [44]: 0.3 mT for Fe 8 and 25 mT for Mni 2 , which 
are close to the measured ones, before these experiments were performed and come (in 
their theory) from the nuclear spins. It was shown that the dipolar fields produce the 
bias ^ H which a few orders of magnitude larger than the value of the tunneling splitting 
and, therefore, can block the tunneling. However, the fast nuclear dynamics (transverse 
relaxation in nuclear subsystem or, in other words, T 2 processes) broadens the resonance 
line and opens a channel for tunneling. This explains the origin of the hole (its width is 
defined by the width of the distribution of the hyperfine fields) in the field-dependent 
relaxation rate (which is proportional to the total distribution of internal fields). 

All these low-T experiments performed on Feg [41,43,45,53,56]] demonstrate the 
square-root relaxation law at short times (about first 100 sec, see e.g. Fig. 17). 
According to the theory of Prokof'ev and Stamp, this law comes from the time- 
dependent distribution of fluctuating internal dipolar fields in a sample. Fast initial 
transitions change the total distribution of internal fields across the sample that can push 
some molecules out of resonance but bring other molecules into the resonance 
everywhere in the sample and allow continuous relaxation. For longer times the 
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experimental relaxation data are better fitted by the stretched exponential law with 
pM).4 [41,43]. 

In Mni2, Thomas and Barbara [57-59,21] show that the relaxation behaves non- 
exponentially even at relatively high temperatures (up to 2.8 K) that can be understood 
as a consequence of intermolecular dipolar coupling. (At low temperatures non- 
exponential relaxation was seen in many experiments. See, for example, 
[4,6,7,11,53,62]). The exponential relaxation can take place only in the limit of non- 
interacting (with each-other and with nuclear thermostat) molecules when we can 
simply write dM(t)/dt=-x " M(t). In the case of interacting molecules, the right-hand side 

of this equation contains some different value M*(t) =f(M(t)) which gives deviation 
from the exponential law (note that a simple distribution of relaxation times gives also a 
non exponential behavior). As it was found in [57], below approximately 1.7 K the 
relaxation follows a square -root law with characteristic time x weakly dependent on 
temperature whereas at the temperatures above 2.5 K the characteristic time follows the 
Arrhenius law, but the relaxation is still non-exponential because of the dipolar 
interaction between molecules. To fit the experimental relaxation curves, it was 
necessary to use the stretched exponential law with (3(T)<1. Below 1.9 K (3(T) is 
approximately constant near 0.5. However, at low temperatures the stretched 
exponential was not really satisfactory. As can be seen from Figure 31, the square-root 
law fits all the experiments up to -1.7 -f 1.8 K (see also Thomas and Barbara in [58]). 

At higher temperatures, between 2 K and 2.8 K, dipolar interactions play a less 
fundamental role, but their influence is quite observable, e.g. from the fact that the 
relaxation is clearly non exponential (the exponent of a stretch exponential fit increases 
from 0.5 to 1, see [57-59,21]) and also by the shift of the maximum of the relaxation 
curve in magnetic field. This shift, found by Thomas and Barbara in Mni 2 at different 
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temperatures [57-59,21], is due to the evolution of the internal field with temperature. 
They suggested [57,21] that this shift can produce a square root like relaxation laws of 
origin different from the Prokof'ev and Stamp law. Instead of being at zero Kelvin, this 
law, with temperature-dependent relaxation times, can hold at high temperature, even if 
the system is equilibrated. This suggests that a quantitative description of the Thermally 
Assisted Regime in Mni 2 requires to take into account interactions with the spins 
environment (not only with the spin-phonon interaction). This also explains [see [21]) 
why the resonance lines have a Lorentzian shape [11,12] although the relaxation is not 
exponential. Note that even if the observed resonances can be well fitted to Lorentzian 
law [11,12] (and not to the Gaussian one associated with hyperfine and dipolar fields 
[60]) this does not mean that hyperfine and dipolar interactions are not relevant. The 
widths of the dipolar fields and hyperfine fields distribution in Mni 2 are comparable, 
and these fields can participate in the resonance broadening but their non-Lorentzian 
character is hidden by the long Lorentzian tails (and experimental error bars). The fact 
that the observed Lorentzian line-shape does not prove that the system of molecules is 
equilibrated, was show in [21] where it is shown that the shape of the resonance 
depends on the history of the sample (field cooled or zero-field cooled, quenched or not 
quenched). Finaly one should mention that the published line- widths depend on the 
authors; they are going from 20 mT [60] to 30 mT [53] or 35mT [12]. In [21] a 
continuous increase of the resonances line-width was observed with the index of the 
resonance at temperatures close to 3 K. This effect being not symmetrical with the 
magnetization state M=0, one cannot infer the dipolar field distribution here. It may 
rather be due to the intrinsic line-width A, which is at the higher temperatures (when the 
tunneling takes place at the top of the barrier) of the order of 0.5 K (which is larger than 
the other contributions [21]) and which increases in average with the applied field. This 
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observation confirms that, even if dipolar interactions and hyperfine interactions (with 
longitudinal Ti relaxation processes, driven by the fast dipolar processes [44,47,61]) 
are essential in the range 2 K to 2.8 K, at the highest temperatures (say between 2.8 and 
3 K) the intrinsic line-width A plays the most important role. In this case i.e. at 
temperatures rather close to the blocking temperature (3.0 K), the transitions are 
homogeneously broadened and the observed Lorentzian line-shapes really attest for 
thermal recovery of the spin-phonons system. (Note, that fast dipolar flip-flop processes 
are also able to provide thermal quasi-equilibrium on each side of the barrier, 
producing, simultaneously, rapidly fluctuating fields acting on each molecule.) In this 
limit the spin-phonon interaction becomes predominant leading to the Thermally 
Assisted Tunneling Regime. 

The observation of history-dependent peak shapes in Mni 2 at relatively low 
temperature (1.6 K) [21], leads us to ask the question of the validity of experimental 
determinations of Lorentzian line-shapes. To extract this distribution from the measured 
relaxation curve vs field, it is necessary to assume an analytical law to fit the data. As an 
example, Friedman et al. in [60] fitted the experimental curve to exp(-t/x). However, as 
this is shown above [57-59, 21,53] the relaxation deviates clearly from an exponential 
even at these high temperatures (2-2.6 K). This must lead to large errors in the 
relaxation time (obtained at each field) and, therefore, in the resonance line-width. Note 
that even in exponential regime, error bars on the equilibrium magnetization M eq in 
M(t)=M eq +(Mi n -M eq )exp(-t/x) lead to uncertainties in the line-shape. The value of M eq 
can be obtained accurately only in the case of very fast relaxation, i.e. either very close 
to the blocking temperature of 3 K (and in this case the observation of a Lorentzian 
form makes no doubt) or in the presence of a large transverse magnetic field which 
increases the tunneling gap [85]. 
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In conclusion, each molecule can tunnel under the effect of fast fluctuating field 
originating from dipolar and hyperfine interactions, and also of spins-phonons 
interactions at higher temperatures. It is only in this last case (where the transition width 
is intrinsic) the spins and phonons is equilibrated (no hole in the spin energy 
distribution). 
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4.2. Thermally Assisted Tunneling Regime. 

The first attempt to apply phonons-based mechanisms was made by Villain et al. in 
[63] and Politi et al. in [64] before the resonant tunneling in Mni 2 was confirmed 
experimentally. Ref. [64] considered the giant spin model for single Mni2 molecule in 
longitudinal magnetic field (Eq. (1) with no S z 4 term) with coupling to the acoustic 
phonons in the form: 

H sp -ph=£ (ri/2N u M u co q )- 1/2 [iqV q (S)C q + - iqV q + (S)C q ], (19) 

where V q (S)~D(S x S z +S z S x ), N u is the number of unit cells, M u is the mass per unit cell 
and C q is the phonon annihilation operator. They have ignored the possibility of the 
tunneling (due to the forth-order anisotropy term) at the top of the barrier (as well as at 
the bottom) and have got the result x _1 ~(SH z ) 3 for the relaxation rate, in contradiction 
with the experiments which demonstrated the existence of a minimum of the relaxation 
time near H z =0 instead of maximum (see Figure 5). Later, this theory was extended by 
including the interactions with nuclear spins (the possibility of tunneling was ignored 
again and only the longitudinal part of hyperfine interaction was considered) [65]. It 
was found that combination of these two mechanisms (strictly speaking, the sum of two 
different curves) can give a minimum of the relaxation time at zero field. This theory 
was also unable to describe the resonant behavior of the relaxation curve. Phonon- 
mediated tunneling relaxation was considered in the theory of Garanin and Chudnovsky 
(see [66]) which involved also random hyperfine fields. However they ignored higher- 
order anisotropy terms in the Hamiltonian and based their calculations on the 
perturbation theory for smaall transverse field H x . All these models are qualitatively or 
quantitatively in contradiction with the experiments. Nevertheless, they are not incorrect 
and they show (directly or indirectly) the importance of the hyperfine interaction as well 
as of the intermolecular dipolar coupling (see, for example, Burin et al. in [67]). The 
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problem is that to get a quantitatively correct answer, it is necessary to consider (even at 
temperatures around 2 K) that the tunneling effect involves all the interactions of the 
spin of each molecule with the environment of spins and phonons. This is very difficult 
task. For example, Dobrovitski and Zvezdin in [68] have concluded that a correct 
description of the jumps width of the hysteresis loop (see Figure 3) requires to take into 
account fluctuating internal fields since they found for pure giant spin Hamiltonian a 
huge discrepancy with the experimental results. However, they just limited themselves 
to the suggestion that the origin of this field could be of dipole-dipole nature. They 
estimated an average value for jumps width using a Gaussian distribution of fluctuating 
field. In some sense similar calculations were made by Gunther in [69] who calculated 
the width of the jump in the hysteresis loop and concluded that it is necessary to involve 
the dynamical transverse magnetic fields to avoid the discrepancies with the experiment. 

The first theory taking into account the fourth order anisotropy terms is of Luis et 
al. [50]. They presented a theory of resonant quantum tunneling of large spins through 
thermally activated states which includes: phonon-mediated transitions between the 
states m and m' with m-m'= 8m=± 1 (in a simplest form S X S Z +S Z S X ), resonant tunneling 
due to fourth-order anisotropy terms, and transverse magnetic field. They assumed that 
the transverse magnetic field originates from the combined action of dipolar and 
hyperfine fields. They concluded that, of course, internal fields can not explain the 
minima on the relaxation curve alone but, together with fourth-order terms, these fields 
can account for the experimentally observed behavior of the relaxation in Mni 2 . To 
obtain the life-time of the excited levels (due to phonon-mediated transitions) they 
applied a standard master equation. Finally, the spin-relaxation rate was averaged over a 
Gaussian distribution of longitudinal dipolar fields (together with hyperfine ones). 
Despite the fact of the inequivalence in treatment of tunneling between the resonance 
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states and phonon-mediated transitions, this model was the first one which described 
(qualitatively) the hysteresis loop of Mni 2 [11], as well as other experimental results at 
high temperatures (T>2.5 K). It is important to note, however, that in this model the 
magnetization relaxes always exponentially "after a brief non-exponential relaxation" 
(see also the same authors in [78]). However, as we have seen above, the magnetization 
relaxes non-exponentially, unless the temperature is very close to the blocking 
temperature (about 3 K). 

Recently Fort et al. in [51] have improved calculations of Villain et al. [63-65] by 
adding tunneling through the top of the barrier. They derived the master equation: 

dN m /dt=£ N m . pYm _ p m +£ N m+pYm+p m -N m £ (Ym m " p + Ym m+p )-(N. m -N m )r m , (20) 

where N m is the number of molecules in spin state lm> , y m p - the spin-phonon relaxation 
rate from state Im) to state lp> and r m is the tunneling relaxation rate from state lm> to 
state l-m) . This equation includes the phonon mediated transitions with 8m=± 1, ±2 
together with tunneling between the states lm> and l-m) . Using Eq. (20) authors 
concentrated on the investigations of the first resonance (H z =0) on relaxation curve (see 
Figure 5). They noted that in the region of the validity of their theory (H z <3 kOe) the 
relaxation rate calculated at the temperatures T=2.8 K and T=2.97 K are in good 
agreement with the experimental one and this agreement between the theory and the 
experiment becomes better at higher temperature. They also noted that a general 
treatment of this problem requires to incorporate other interactions (dipolar, random 
fields, etc.) since the model with fourth-order anisotropy term gives only transitions 
with 8m= ± 4. To modify the theory, it was suggested to include some transverse fields 
(of any nature) which obey selection rule 8m=± 1 (i.e., contain terms like S x or S y ). 
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Very recently Leuenberger and Loss [52] presented a theory of the relaxation in 
Mni2 at high temperatures (2 K and higher) based on thermally assisted spin tunneling 
in weak transverse magnetic field. They solved the standard master equation (for the 
reduced density matrix p(t)) which includes both the resonance tunneling due to fourth- 
order anisotropy term (as well as transverse fields) and phonon-induced transitions with 
8m=± 1, ±2. The origin of the transverse field was attributed to a misalignment of 
m =l° between the field direction and easy axis, extracted from the experiment of 
Friedman et al. [60] (see, also [4]). As it is noted by authors, this model is in 
"reasonably good agreement with all the experimental parameter values known so far". 
The differences with [51] are the following: 1) a more general spin-phonon interaction 
was considered; 2) transitions due to a transverse magnetic field were included; 3) 
longitudinal fields were not limited to the resonance near H z =0. All that allowed them to 
get an independent description for each resonance in the experimental curve of Figure 5. 
To obtain a continuous description of the relaxation time vs H z the authors applied the 
Kirchhoff's rules by associating with each independent path from 1-10) to 110) the 

probability current J n =dp n /dt, where p n is the reduced density matrix p(t) for particular 
path n. This theory can effectively be applied to Mni 2 , but after some correction given 
below. [52] allows to recover the observed Lorentzian line-shape for the resonance 
peaks, but the procedure that was used to get the experimentally observed line-width 
was "confusing". The authors cut off the calculated Lorentzians to the appropriate 
values since these last (the Lorentzians) were found extremely high and narrow. This 
procedure, of course, was supported by mathematical arguments, but the problem is that 
the authors have missed (as well as Fort et al. in [51]) in their theory most important 
contribution, namely, the value of the tunneling splitting in denominator of their 
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formula for the tunneling rate. As we will discuss below, exactly this missed term gives 
most important contribution to the width of the resonance line. 

It was also shown in [52], that even and odd resonances should have different size 
since even resonances are induced by S+ 4 or S. 4 , whereas odd resonances are induced by 
combinations like S X S+ 4 or S X S. 4 (similar effect was found experimentally by Thomas 
and Barbara [11,21]). In addition, based on experiments of Caneschi et al. in dilute 
samples [70] (compared with the same measurements on powder sample), authors 
completely ignored any dipolar fields as well as any hyperfine couplings since its 
contribution should give rise to the Gaussian distribution, whereas the peaks of 
relaxation rate, observed experimentally, have the Lorentzian shape. However, we have 
noted already (see the next to last paragraph before Section 4.2.), that it is not so easy to 
conclude from the present experiments (at T<T B ~3 K) about the real shape of these 
peaks (even if they are the Lorentzians) close to tails. Even if the phonons dominate in 
the relaxation mechanism at T>2 K (with no doubts), the influence of dipolar and 
hyperfine couplings can not be negligible because (at least) the relaxation behaves non- 
exponentially at these temperatures (see [21,53,57-59,62] and references therein). 

In the following we would like to give a qualitative picture of the Thermally 
Assisted Tunneling phenomenon including couplings to the environment (with respect 
to the resonance line- width). We do not state that our calculations are complete and 
want just to show how this works. (Note, that the similar calculations have been done by 
Prokof ev and Stamp several years ago but they did not publish it) 

The Hamiltonian of interest for the Mni 2 molecule can be written as follows: 

H=H G + H sp .ph + H h yp+ H dip , (21) 
where Hg is given by Eq. (1) and we take spin-phonon Hamiltonian in the simplest form 
of Eq. (19) with V(q)=D(S x S z +S z S x )=(D/2)[(S + +S.)S z +S z (S + +S.)]. The third term in Eq. 
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(21) describes the hyperfine interaction of central spin S with nuclear spins Gk of each 
Mn ion (i.e., N=12 nuclear spins): 

H hyp =l/S£ (Tit0k/2)So k . (22) 

k=i 

And, finally, the last term gives the dipolar intermolecular interaction: 

H dip =l/2£ VijSiSj, (23) 

where Vjj depends on cube of the inverse distance between molecules. Let us assume 
now that the system is close to the characteristic magnetic field H n (Eq. (3)) when some 
pairs of levels from the two opposite sides of the barrier come into resonance. As this 
pair of levels is well separated from the others, we can truncate the giant spin 
Hamiltonian H G to a two-level one H m>n . m and consider the pair of levels lm)and ta- 
rn) with the corresponding energies of E m ° and E n . m ° (E ra °= -Dm 2 -K,|m 4 ): 

Hm,n-m = Am,n-m%t ^m,n-m^Z5 (24) 

where A m , n _ m is the tunneling splitting (due to fourth-order anisotropy terms, transverse 
components of dipolar and hyperfine interactions and transverse components of external 
magnetic field which can be originated in the experimental misalignment of the crystal), 
^m,n-m is longitudinal bias: 

^ m , n . m =(2m-n)[(E m °-E n _ m °)/(2m-n)+l/sf; ( ^ oy2)a k z + £ (V 0d /2)mj -gjj. B HJ/2 (25) 

*=1 7*0 

and x x , x z are the Pauli matrixes. In Eq. (25) the third term comes from the longitudinal 
part of dipolar interaction, where m/ is the spin state of j-th molecule. In this case we 
can get for the maximum value of the tunneling probability with no phonons: 

Pm,n-m = ^m,n-m /(^m,n-m + ^m,n-m ) (26) 

The longitudinal hyperfine couplings give a Gaussian spread to each giant spin energy 
level E m . This means that each energy level E m is actually split into a Gaussian multiplet 
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with N+l different polarization groups of N=12 nuclear spins. In reality one should 
consider tunneling between sublevels from opposite sides of the barrier, which are in the 
resonance for a given field. However, in the simplest case (zero approximation) we will 
consider the longitudinal internal bias (which comes from the hyperfine coupling) just 
as some "mute" variable 8 with the Gaussian distribution: 

G(E)=(27C0 2 )" 1/2 exp(-8 2 /2o 2 ). (27) 
In this case the half-width o ~N 1/2 (Oo with Cflo=( 0\) and according to the experimental 

results given above [53,21,87] o ~6 mT (i.e., O)o = 1.75 mT). As for the dipolar 
coupling, we do not know it reliably (except from some preliminary measurements 
which give E D =2o o ~20 mT) and, therefore, for the distribution of dipolar fields we take 
the same value o ~6 mT. A remark should be made here. The distribution of dipolar 
fields strongly depends on the value of the magnetization (e.g. M~0 after zero-field 
cooling or M ~ M s after field cooling, where M s is the saturated magnetization) and on 
the shape of the sample. Only at sufficiently strong annealing M < 0.5 IM S I one can 
observe the Gaussian distribution for dipolar fields. In the following we assume zero- 
field cooling (an annealed sample). Thus, we can rewrite Eq. (25) in a more transparent 
form (e includes the contributions from both hyperfine and dipolar fields): 

^ m , n . m =(2m-n)[g^i B (H n -H z ) +e]/2, (28) 
where H n is the characteristic field from Eq. (3) (actually H n are calculated using Eq. 
(1)). Choosing a Gaussian distributions for both dipolar and hyperfine fields, and with 
the same o , is a great simplification which leads to think that dipolar fields are more or 
less ignored, since they behave as "non-interacting" hyperfine fields. Strictly speaking, 
we are not far from this suggestion in our simple model, but here we do not see the 
reasons to play with the form of distributions. In reality, however, one should include 
also the flip-flop transitions between molecules which come from the transverse part of 
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dipolar interaction. Moreover, when system relaxes, the total bias field (internal plus 
external) changes in time. The flipping of the molecules produces time-dependent 
fluctuations of dipolar fields which also cause the transitions in the nuclear subsystem. 
Therefore, each molecule feels rapidly fluctuating field ^=^(t) [44,47,61,72,73]. This 
means that the Boltzmann distribution is time-dependent and an additional time- 
dependence produces deviations from the exponential law (and these deviations are 
increasing with decreasing of the temperature). Ignoring all that, we assume here that all 
these processes are fast enough to keep the thermal quasi-equilibrium at each side of the 
barrier with the "static" Boltzmann distribution. This leads, of course, to the exponential 
relaxation, but the above simplification is sufficient for our purpose which consists in 
estimating the width of the resonance line in the case of annealed sample. 

The transverse part of the hyperfine interactions contributes to the Berry phase of a 
central spin in producing a random complex phase (for a details see Tupitsyn et al. 
[23,31] and Prokof'ev and Stamp [47,72,73]). To our present concern, one can say that 
the transverse hyperfine interaction will act as a transverse field, and therefore change 
the tunnel splitting to an effective one: 

A m , n . m =A m , n . m ($), (29) 
where A m>n . m (0) is the tunneling splitting in zero external transverse magnetic field with 
no hyperfine interactions and $ is an additional phase from hyperfine interactions. It 
would be wrong to say that the distributions of internal transverse and longitudinal 
fields are always the same. However, in zero approximation (and for strongly annealed 
sample) we can use this suggestion and consider the Gaussian distributions with the 
same width Oo for both transverse and longitudinal fields (of course, it is easy to take 
different Go for these fields). 
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Next, we should take into account the magneto-acoustic interaction. It was shown 
by Kagan and Maksimov [74] that the correct contribution (to all orders in A) from 
inelastic phonon processes to the transition rate can be written as follows: 

rm,n-m = 2A m n . m W^ m /(^ mjn . m + A m n . m +T) W m ), (30) 

where r|W m is the phonon line-broadening. Note that this equation clearly shows the 
expected Lorentzian form vs a longitudinal fields for a system at quasi-equilibrium (see 
Eq. (28)). The most important transition here is the transition to the next upper level 
(from Im) to lm-1) ) with phonon absorption. From Eq. (19) with V(q)=D(S x S z +S z S x ) 
we get (see also [63-65,90]): 

W m (1) =(3D 2 (S+m)(S-m+l)(2m-l) 2 E m . 1 , m 3 )(87ipc 5 Ti 4 [exp(E m _ 1 , m /T)-l])" 1 , (31) 
where E m _i >m =E m _i-E m , p=M/a 3 is the mass density (a is the lattice constant), 
c=(kB/r|)®D(Vo/67l 2 ) 1/3 (see [75]) is the sound velocity with the Debye temperature D 
and V is the unit cell volume. According to recent measurements of specific heat in 
Mn 12 by Gomes et al. [76] D =(38±4) K and from [1], we can take V =3716 A 3 and 
p=1.83xl0 3 kg/m 3 . 

To see the dominant contribution to the Lorentzian line-width (see Eq. (30)), one 
must compare £ m ,n-m, A mjn . m and r|W m . As already noted, the estimated width of the 
Gaussian distribution of internal fields is approximately 2ao=12 mT (for hyperfine 
interaction, at least). Let us concentrate for the moment to the high temperature regime, 
with thermally activated tunneling from levels m smaller or equal to 4. In zero field A m> . 
m (0) gives l.lxlO" 2 K and 0.34 K, for m=4 and m=2 respectively, whereas r\W^ l) 
« 1.02xl0" 5 K and r|W 2 (1) =6.08xl0~ 7 K (for T=2.6 K and c=1.4xl0 3 m/s). It is only in 
the region of m=6 or larger the values of rjW m (1) and A raj . m (0) become more or less 
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comparable. In order to take into account odd values of m (as well as nonzero n), we 
have to include internal transverse fields (as it is discussed above Eq. (29)). 

We can also consider the phonon-assisted transitions with 8m=± 2. In this case we 
take for simplicity V(q)=D(S x 2 - S y 2 )= D(S + 2 + S. 2 )/2. For transition from lm> to lm-2) , 
this yields (see also [52]): 

W m (2 k3D 2 (S+m)(S+m-l)(S-m+2)(S-m+l)E ra . 2 , m 3 )(87ipc 5 i 1 4 [exp(E m . 2 , m /T)-l])- 1 . (32) 
This equation gives r|W 4 (2) = 2.99xl0" 5 K and r|W 2 (2) = 1.16xl0" 5 K (also for T=2.6 K). 
These numbers show that transitions with 8m= ± 2 can not change the situation with 
relative contributions of the values r|W n and A mj . m (0) to the Lorentzian line- width, 
however, we will use the average value W m =(W m (1) +W ra (2) )/2 in Eq. (30). Of course, 
spin-phonon interaction contains other terms but it is unlikely that they will give 
contribution which is of the order of magnitude larger than W m . 

This means that in thermally activated regime (near the top of the barrier) the line- 
width is defined mainly by the tunneling splitting of the resonant levels and by the 
internal longitudinal fields (see Eq. (28)), but not by the phonon line-broadening. Note, 
however, that the phonons play an essential role in "linking" the states on the same side 
of the barrier (otherwise only the ground-state would be occupied, unless one admits 
that all the bias is dynamical). The dipolar flip-flop processes also can cause transitions 
between the energy levels providing thermal equilibrium. 

We should note that Villain and Fort et al. in [51] as well as Leuenberger et al. in 
[52] have slightly different definitions for the tunneling rate, but which in general, look 
like: 

r m , k =2A m , k 2 W m /(^ m , k 2 + r| 2 W m 2 ), (33) 
where m and k are the levels in the resonance and ^ m>k is the bias of Eq. (25) (or Eq. 
(28)) with no internal bias fields. According to this formula, the half-width of the 
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resonance peak should be of the order of r|W m which is actually very small in 
comparison with A m ,. m in zero field for upper levels where the tunneling takes place at 
T>2 K. (see above, just after Eq. (32)). The reason of this problem, is that these authors 
missed the term relative to the intrinsic width A in their tunneling rate (see Eq. (30) and 
(33)) which is absolutely not negligible for thermally activated tunneling. 

In the thermally activated regime, the value of the tunneling rate can be evaluated 
for each resonance (i.e. each value of n) by taking the product of the intra-well 
(Boltzmann) and inter- wells (tunneling) transition probabilities. This has to be summed 
up to all the contributions from different lm> : 

x- 1 n (H z )=r 1 (H z ) £ r m , n . m exp[(-E m °-g^i B H z m) /k B T], (34) 

m 

where Z(H Z ) is the partition function. (For simplicity we omit here the time dependence 
of the Boltzmann factor, i.e. we still assume quasi-equilibrium). The plot of x _1 n (H z ) 
gives the expected Lorentzian line- shape of resonance peaks. The width is rather 
sensitive to parameters such as the distribution of internal fields or the sound velocity. 
However realistic values for these parameters allow to obtain the experimentally 
measured line-widths. To give some numbers, we obtained the width of about 25 mT 
(see [60]) with o = 6 mT at T=2.6 K. To get the same height, it was necessary to 
change only one parameter from all the set. We put the sound velocity to the value of 
c=1.4xl0 3 m/s, the other parameters are the same as given above (see Eq. (1) and Eq. 
(31)). Our curves show slightly faster decay near the tails (which is the consequence of 
the averaging over a Gaussian distribution of fluctuating internal fields). However, if we 
neglected these internal fields, Eq. (34) would give too narrow and sharp resonance 
lines. 
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Before concluding, a few remarks should be done. We limited ourselves by taking 
an average of Eq. (30) over fluctuating internal fields (longitudinal in bias £ m , n -m and 
transverse in A m , n _ m ) but this is just zero approximation, to show that even this 
approximation can lead to correct resonance line-width. As we noted already, in reality 
the interaction with the nuclear subsystem spreads each giant spin energy level into a 
Gaussian multiplet and one should consider all possible transitions between the resonant 
levels from the opposite sides of the barrier inside such multiplets. Moreover, the 
internal bias field e(t) is actually time-dependent (i.e., it varies not only in space) and 
produces, in general, a time-dependent Boltzmann factor (Z" 1 (H t (t))«(exp[(- 
E m °+mH t (t))/k B T] where H t (t)=g(i B (H z -£ (t)) and, as a consequence, deviation from 
exponential relaxation occurs. 

To conclude our discussion of the high-T behavior of the relaxation in Mni 2 - 
acetate, we should emphasize the following main results: 1) At relatively high 
temperatures (T>2 K) the relaxation is dominated by the Phonon- Assisted Tunneling 
Mechanism; 2) The Lorentzian shape of the resonance peaks is determined by Eq. (30), 
describing the tunneling between the levels which are in resonance in a given magnetic 
field; 3) Despite the dominant role of phonons the resonance line-width is determined 
mainly by the value of the tunneling splitting A m n . m together with internal bias 8 (Eq. 
(28)) which originates from the hyperfine and dipolar interactions; 4) All mentioned 
theories predict exponential relaxation of the magnetization while the experiments show 
a non-exponential behavior at these temperatures (see [4,6,7,57-59,21,53,62]). This fact 
clearly shows the important role played by intermolecular interactions as well as by 
interactions with nuclear spins which both produce time-dependent bias field causing 
the non-exponential behavior of relaxation. 
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4.3. Ground-state Tunneling. 

Let us now concentrate on the low-temperature limit of T«T C where tunneling 
takes place on the ground-state only. In the instanton approach (see below), we assume 
that k B T«Q , where Q is the "bounce frequency" of the instanton transition (which is 
roughly the distance from the ground-state level to the first excited one. In this case only 
the lowest levels are thermally populated and we can truncate the "giant spin" 
Hamiltonian (Eq. (1) or Eq. (12)) to an effective one (which is valid only at energies « 
Q. ), describing two states IS Z )=±S and their mixture, separated by the tunneling 

splitting A«Qo- Moreover, we should couple the "giant spin" S of each molecule to the 
nuclear "spin-bath" {a k } with k=l,2,...,N. In the case of hyperfine coupling, the latter 
can be described by Eq. (22) with u) k «Qo ((Qk»A in most cases) and with N equal to 
actual number of nuclear spins inside molecule (one can take N=12 for Mni 2 neglecting 
the effect of hydrogen. . ., and N=146 for Fe 8 , including 120 hydrogen, 18 nitrogen and 8 
bromine atoms, neglecting the few percents of Fe57). Without hyperfine interactions, 
the total nuclear spectrum, containing 2 N states, is almost completely degenerated, with 
only a tiny spreading ~T 2 A of levels caused by the inter-nuclear dipolar interactions. 
With the hyperfine interaction, the nuclear levels spread into a Gaussian multiplet of 
N+l polarization groups around each giant spin level (see [47,61,72,44,23]). The half- 
width of this distribution a is of the order of (OqN 112 (where COo=(l/N)^^ (Ok). Note 

that, in each polarization group, the hyperfine levels are also distributed according to a 
Gaussian half-width ~T 2 _1 . Typically the different polarization groups completely 
overlap within of the Gaussian envelope of Eq. (27). The nuclear T 2 -processes 
(transverse relaxation) are responsible for transitions inside each polarization group, 
whereas the nuclear Ti -processes (longitudinal relaxation) provide transitions between 
different polarization groups. According to [45], the half- width of the hyperfine 
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distribution in Fe 8 is of the order of 0.6 mT, yielding cOo~0.05 mT (for the estimate we 
take only hydrogen atoms with a=l/2, i.e., N=120). Since cOo»T 2 " 1 (typically T 2 " 1 ~10~ 7 
K), the nuclear spin dynamics is slaved by dynamics of S. We also suggest that the TV 
processes are long and therefore not relevant, at these temperatures. The reason is that 
the Ti-processes are driven by the dipolar flip-flop transitions, which are essentially 
frozen at T«T C . In general, we should consider two effects: (i) the effect of nuclear 
spins on the giant spin dynamics during the tunneling and (ii) the effect of the motion of 
S on the nuclear spins. Both effects have to be handled in a self-consistent way. 

Using the instanton technique, we can write the effective Hamiltonian for a single 
molecule in the following way [23]: 

N 

H eff (1) =[2A x.cos(tiS - (3 n «H + ^ o^n *tf k ) + H.c] 

k 

N _ N 

+(1/2)[t z £ Cflk"l k «o k +£ a^mk'Gk] 

k k 

+t V k i ap a k a a, p , (35) 

k*l 

where x describes the "giant spin" of molecule (x and a are both the Pauli matrices). 
The first term in Eq. (35) is non-diagonal term (due to x+) which operates during 

transition of S. It produces a time-dependent field 7k=(a)kS/2S) acting on each a k , and 
this leads to a k to flip. If we expand out the cosines, we see that we have a whole series 
of terms like ~T±r a ?,y8t3ki a <5k2^<3k3 y <5k4 '• • • in which an instanton flip of the giant spin 
couples to many different nuclear spins simultaneously, i.e., a single instanton can 
simulate multiple transitions in the nuclear bath. The probability that a k will flip during 
a single instanton passage between two quasiclassical minima ISi>and IS 2 )is loc k l 2 /2. 
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Thus, the average number of nuclear spins that will flip each time when S flips (so- 
called "co-flipping" amplitude) is approximately: 

Ml/2)£ Icy 2 (36) 

k 

which, in principle, can be larger than 1. One can easily calculate the dimensionless 
constants 0Ck and (3o in the case of the simple bi-axial Hamiltonian of Eq. (6) with H hyp 
given by Eq. (22). The answer is [23,31]: 

oc k n «a k » (TiCDk^Qo) [-io y +(D/(E+D)) 1/2 o x ] (37) 

pen -H « (7lg^ B S/ao)[-iH y +(D/(E+D)) 1/2 H x ] , (38) 

where n is a unit vector in the (xy) plane and the "bounce frequency" Q.q reads as: 

Q = 2S(ED) 1/2 . (39) 
(For non-zero values of K|_ see Eq. (14-18)). Eq. (37) tells us that A<1 in the particular 
case of E and D as in Fe 8 (or Mni 2 , note that non-diagonal fourth-order terms 
renormalize Eq. (37,38) but for Ki./k B =-3.28xlO" 5 K it gives a k ~((O k /Q. C 1 ) with 
Cl«1.56). 

The second term in Eq. (35) is diagonal, which operates when S is in one from its 
two qusiclassical minima. Let us introduce two corresponding fields y k (1) and y k (2) . In 
general case (any non-zero external and internal magnetic field) Si and S2 are not 
antiparallel. It is easy to see that the sum and the difference between these two vectors 
define uV' and 0\ L i.e., 

0)k"lk=Yk (1) -Yk (2) (40) 
co^m^Yk^+Yk^, (41) 
where l k and m k are mutually perpendicular unit vectors. The longitudinal coupling oV' 
gives the change in energy of a k before and after S flips (i.e., the difference in effective 
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field acting on G k before and after transition of S). The transverse coupling Cflk defines 
the deviation of initial and final orientations of S from the ± z-direction (which is easy- 
axis for the Hamiltonian of Eq. (12)). In the case of biaxial Hamiltonian for small values 
of H x (g(i B H x «S(E+D)) we get: 

(Ok" lk'tfk ~®k*Oz (42) 

©k 1 m k -Gk~(aW B H x /(2S(E+D))) -o x (43) 

Finally, the third term in Eq. (35) describes very weak internuclear dipolar 
coupling (IVki <x ' 3 l~T2" 1 ). To complete our Hamiltonian, we should also include the 
dipolar-dipolar interactions between molecules. Ignoring transverse part of this 
interaction (which leads to flip-flop processes), we get: 

H D =(l/2)£ V m , v (D) t/-t z v , (44) 

where IV^ V (D) I~1 mK. Thus, to work in the low-T Quantum Regime, one may use the 
following effective Hamiltonian: 

Hef t -H eff (1) +H D . (45) 
We do not include the spin-phonon interaction into this Hamiltonian since the phonons 
can play no role at T<T C . 

Now, after the effective low-T Hamiltonian is established, we would like to 
discuss briefly some effects which follow from the coupling of the central spin to the 
spin bath (for detailed explanation one should read original papers [72,47] and review 
[73,61]). 

Let us start from the nondiagonal term in Eq. (35) (the first term). Since 0Ck is 
complex quantity, one can expect two effects. The imaginary part of 0Ck gives a 
renormalization of the effective tunneling splitting depending on the coupling constants 
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(Ok. One can expect an increase (decrease) of the effective tunneling splitting with u)k. 
The width of the distribution of hyperfine fields also depends on C0k. It should also 
increase with (Ok. The real part of 0Ck adds an extra random phase to the Haldane phase 
Pon»H. As we have seen already, the average number of nuclear spins that flip together 
with S is proportional to loCkl 2 . These flips modify the total phase of the bath state and, 
consequently, randomize the phase of the giant spin (between the instanton and anti- 
instanton) producing phase decoherence of the tunneling process that can completely 
block the latter (this is called, "topological decoherence'"). Indeed, we know that a half- 
integer spin can not tunnel in zero transverse field. We can imagine that in different 
molecules the transition of S is accompanied by the different number of the nuclear 
spins. If the total flipping spin (S together with the nuclear spins G k ) is integer, the 
tunneling is allowed. Otherwise, the tunneling is blocked. Thus, the possibility to tunnel 
depends on the particular environment state inside of each molecule (nuclear spins), and 
over the entire sample the transition of S can happen at random. 

The second term in Eq. (35) produces an internal bias field (we put here (Qk 1= 0) 
£=(1/2) ^ ©k"<7k Z acting on S. Together with the first term it gives the Hamiltonian of 

the biased two-level system (see Eq. (24)) with an effective tunneling splitting 
2A$=2A cos($) where 4> is the complex phase (we assume here that there is no nuclear 
spin dynamics itself, i.e., I Vki a ' 3 l=0). In this case the tunneling probability is given by: 

P (0) (t)=(4Ad> 2 /E 2 )sin 2 (Et), (46) 
E=±(e 2 + 4A$ 2 ) 1/2 . (47) 
Since £»A_i ,io (in our case Cflo»A_io,ioX only a small fraction of molecules is not 
pushed away from the resonance by the additional longitudinal field £ (£ depends, of 
course, on the particular environment state) and, therefore, are able to tunnel. To 
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estimate this small number of molecules that are close to the resonance, we should take 
the average of P ( } (t) over the ensemble of the molecules with different £ weighted by 
the Gaussian distribution of Eq. (27). This gives (in what follows, for nuclear spins we 
use notation ^ instead of o ) for ^ »Ao (see Chapter 4 in [72]): 

P (0) (t)~p£ W4A*t), (48) 

which oscillates as the Bessel function J2k+i(z) with the amplitude equal to: 

p=(27l) 1/2 Ad>/^>. (49) 
The value of p estimates the fraction of molecules that are able to tunnel. It is easy to 
see that p«l. This mechanism (called "degeneracy blocking") very effectively can 
block the tunneling. However, if we include the dynamics of the nuclear bath 
(IV k i a|3 l ^ 0) , the situation changes dramatically. Due to the interaction between the 
nuclear spins, bias energy 8 becomes time-dependent: within each polarization group 
£(t)=£+8e(t) passes over all the energy range ~T 2 ' 1 and this gives to the system a 
resonance window. Inside this window the total bias field fulfills the condition 
^=^H+e(t)<A (where ^ H is external bias field) and central spin can tunnel. (The same 
mechanism can also destroy coherence by pushing molecules away from the resonance 
window.) Note that tunneling can take place only between the polarization states M and 
-M (M=S Z ) because the energy of the final state E f should be in the resonance with the 
energy of the initial state E ; (lEf-Ejl should not exceed A , at least). This means that if 
initial polarization state of molecule is M, 2M nuclear spins flip when S. As we 
discussed above, the average number of nuclear spins that will flip together with S is 
X«l (for Feg and Mni 2 ). In fact, better possibility to tunnel is realized for molecules 
with initial polarization state M=0. Molecules with M in ^0 also can tunnel (with 
8M=2M) but the contribution into the statistics from these events falls very rapidly with 
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the increase of M. Actually A*(M) falls as ~(X M /M!) 1/2 for M»X (see Chapter 4 in 
[47]). 

If (Ok 1 ^ 0, there is a transverse magnetic field acting on environmental spins. Due 
to this field the initial and final directions of the nuclear spins are not parallel (or 
antiparallel) to each other. Assume that all nuclear spins are initially aligned in 7^ 
(before S flips). After S flips, nuclear spins Gk are not parallel (antiparallel) to the new 
field y 2) acting on them (note that the instanton flip of S is a sudden perturbation for 
nuclear spins, who undergo a non-adiabatic transition). This new state is not an 
eigenstate of the Hamiltonian and nuclear spins must relax in making transitions to 
avoid misalignment with 7^ (this transitions transform their wave-function to the exact 
eigenstate). Thus, the tunneling of S can be suppressed (depending on how slow are the 
transitions in the particular environment state) since the initial and the final states of the 
nuclear bath are not exactly orthogonal. This mechanism is known as "orthogonality 
blocking" mechanism (see Chapter 4 in [72,47]). 

All these effects can be handled by three different kinds of the averaging 
procedure and the final answer can be obtained by combining them depending on its 
importance in each particular case (for analytical expressions see [72,47,73,61]). We 
have seen also that dynamics of the nuclear bath is extremely important since, 
producing rapidly fluctuating hyperfine field, it can help the system to find a resonance 
window by "scanning" over all the range of the bias energy (which for single molecule 
is of the order of T 2 _1 ). For the particular case of Fe 8 or Mni 2 we have X«l which 
means that, in general, molecules in the resonance window relax incoherently (only 
molecules in nuclear polarization state with M=0 can relax coherently in this case). The 
relaxation rate for such incoherent process is given by [47]: 

x N 1 (^)«To- 1 exp(-l^l/^o), (50) 
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To" 1 -2A. 10 ,io 2 /7C 1/2 ^o. (51) 
However, these equations describe just the initial stage of the relaxation. Continuous 
relaxation requires to bring more and more molecules in resonance that can be provided 
by fluctuations of the dipolar field across the sample [44]. When the spin of a given 
molecule flips, it produces the time-dependent magnetic field (long-range) which can 
push some molecules away from the resonance and bring other ones into the resonance 
window depending on the nuclear bath state of each molecule. 

In order to investigate the problem of the relaxation in quantum regime (ground- 
state tunneling), Prokof'ev and Stamp in [44] introduced a kinetic equation for 
distribution function P a (£,r,t) which gives the probability to find a molecule at position 
r with polarization a=±l (i.e., IS Z > =±S) having a bias energy t, at time t. This equation 
reads as: 

dP a (£,r)/dt=-T N '©[Pa&r)- P. a £,r)] - £ J" dr'Qo" 1 J" d£ x N \%) 

a' 

x[P aa . (2) (^,^,r,r')-P aa . (2) (^-aa'V (D) (r-r'),^',r,r')], (52) 
where P aa ' (2) (^,^',r,r') is a two-molecule distribution which gives the probability to find 
a second molecule with polarization a' and bias if the first molecule is with a and 
bias ^, correspondingly. The quantity Q is the volume of unit molecular cell, V (D) (r) is 
the longitudinal part of the dipole-dipole interaction (see Eq. (44)), and integration Jdr' 

is performed over the sample volume. The first term of this equation describes the local 
tunneling relaxation whereas the second (which is analogous to a collision integral) 
describes the influence of the dipolar field produced by the spin flip of a molecule at site 
r' . The solution of Eq. (52) (analytical or numerical) gives the magnetization M(t) as a 
function of time by the following obvious equation: 

M(t)=j dgj drW^P^iO-P.&r)]. (53) 
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If at t=0 the sample is fully polarized, at short times the solution of Eq. (52) can be 
found analytically (for ellipsoidal shape, at least). At the beginning of the relaxation the 
number of the flipped molecules is small (M(t)/M s «l, where M s is the saturated 
magnetization) and, according to Anderson (see [55]), the field distribution of the 
randomly placed dilute static dipoles can be described by the Lorentzian: 

P„(5)=[(l+ccM(t))/2] • [(r d (t)/7i)/{(^-aE(t)) 2 +r d 2 (t)}]; (54) 

r d (t)=(47l 2 E D /3 5/2 )[l-M(t)]; (55) 
E(t)=cV D (l-M(t)), (56) 
where c is the sample shape-dependent constant and V D is the strength of the dipolar 
interaction (V (D) (r)=V D Q [l-3cos(e)]/r 3 ). In this case (short times, the Lorentzian 
distribution) the two-molecule distribution function becomes factorizable (i.e., 
P (2) (1,2)=P(1)P(2)) and Eq. (52) gives the square-root relaxation law [44]: 

M(t)/M s =l-(Tsh 0rt 1 t) 1/2 , (57) 
where Xshort" 1 is the relaxation rate: 

W^riA^-io 2 P(£d)/ti (58) 
with a sample shape-dependent constant T] and the normalized distribution of dipolar 
fields in a sample is P(^d)- When the number of the flipped molecules becomes large 
enough (larger than 10 15 %), the field distribution of the flipped spins becomes non- 
Lorentzian and Eq. (52) should be solved numerically since Eqs. (54-56) are no longer 
valid. This equation also can be solved numerically for nonsaturated sample (i.e., 
Mi n /M s <l, where M in - initial magnetization) as well as for samples of different 
geometry with M in =M s . The latter was done in [44] by kinetic Monte Carlo simulations. 
The main result of these simulations is that the short-time relaxation still obeys the 
square-root law with a sample geometry-dependent constant r\ (very recently this result 
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was confirmed by numerical calculations of Cuccoli et al. in [81]). As for the relaxation 
in nonsaturated sample, very interesting result was obtained experimentally by 
Wernsdorfer et al. in [45] in Fes and by Chiorescu et al. in Mni 2 [85]. They found that 
the short-time square-root law is accurately obeyed for both saturated and non- saturated 
samples. Moreover, the square-root law for a strongly annealed sample or in the 
presence of strong transverse field is even more pronounced, as it can be seen from the 
comparison of the Figure 18 (M in =M s ) with Figure 32 (M in =0)! The theory of Prokofev 
and Stamp (as it is) does not predict the short-time square -root relaxation law for 
annealed samples. Among other possible suggestions, one could say that the distribution 
function P aoc > (2) (^,^',r,r') is also factorizable if M in is small, and all the difference with 
the case of the saturated sample is in numerical constant r|. Anyway, the square-root law 
for annealed sample was confirmed in [25] by the comparison of the tunneling rates 
extracted from the relaxation experiments using Eqs. (57-58) and from the experiments 
based on Landau-Zener model [82,68,69]. These experimental results require further 
theoretical investigation. One should nevertheless mention that this square root 
relaxation is necessarily a short-time regime and a cross-over to another relaxation 
regime must be observed at long times or/and high temperatures. Such a cross-over has 
been observed in Mni 2 for both M->M S [57] and M->0 [85]. The new regime is 
exponential and corresponds to phonons recovery. 

As we understand now, at low T the short-time square-root relaxation law is 
explained by dynamic dipolar interactions (but, to find a resonance window at the 
beginning of the relaxation, the system needs dynamic hyperfine interactions). The 
influence of dipolar fields can also be investigated by measurements of the tunnel 
splitting A(H X ) versus a transverse magnetic field. In Fe 8 at low T, this quantity depends 
strongly (near the nodes) on the value of the initial magnetization (as it was found 
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experimentally in [83] using Landau-Zener method), i.e. it depends strongly on the 
strength of dipolar interactions between molecules (see Figure 33). In what it follows 
we would like to show how to analyze this dependence using very simple language [84]. 

The tunneling probability in the Landau-Zener model depends on the sweeping 
rate of the longitudinal field H z in the following way (see, for example, [82]): 

Pm,m'= 1 -exp(7tA ra , m < 2 (^||,^)/ T|D), (59) 

1/2 

where t>=g|i B (mm') dH z /dt, dH z /dt is the constant sweeping rate and, A mjm - depends on 
dipolar bias fields (^n,^±). At large sweeping rates P miin - ~7tA m;m 2 /r|D. Thus, the average 
probability over the distribution of the dipolar fields reads as: 

Pm,m'(toM7l/TlD)j d£| J d^G(^„,^t ) IA m , m <^,,^)| 2 , (60) 

where G(£n,£i to) is the distribution of the dipolar fields in a sample. We have calculated 
numerically G(^n,^j_t ) for different sample geometries (sphere, cube, parallelepiped) 
and for different initial magnetizations (M(to)=M in ). In the case M in «M s (e.g. zero-field 
cooled sample) the distributions are the Gaussians along all the axes (x,y,z). Then, for 
dH z /dt»l Eq. (60) can be simplified to (for m=S and m'=-S): 

Ps,-s(to)«(3i/T,D)G||(H z )J d^ x J d^ y IA s> - S £iFlU x £ y )l 2 G^ x £ y ); (61) 

G ± (k£ y ) « (27lE DX E DY )" 1 exp[- { (^ x -H x ) 2 /(2E DX 2 )+(^ y -H y ) 2 /(2E DY 2 ) } ] ; (62) 
G l |(H z )-(27cE DZ 2 )- 1/2 exp[H z 2 /2E DZ 2 ], (63) 
Note, that a molecule can tunnel only in the total bias field ^i|-H z <^o/g|i B S. Using this 
condition we put (approximately) £n~H z . Next we would like to calculate Ps,-s(to) near 
the nodes of A s ,_s(H x ,H y ). Let us call y the integral in Eq. (61). Eq. (35) gives the 
effective tunnel splitting: 

A s ,-s=A cosh[7cH y /T y +i7cH x /T x ] , (64) 
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where T x is the oscillation period along the x-axis and T y is the "related" period along 
the y-axis. For Eq. (12) we get: 

T x /T y =[(D+E)/D] 1/2 . (65) 
Near the node (the first node, for definiteness) we can put: 

H X =H X °+8H X +H x dip , H y =H y °+8H y +H y dip , (66) 
where (H x °,H y °) is the position of the node in the (H x ,H y ) plane, 8H x y is the distance 
from the node and H XiY dip is the contribution from the dipolar fields (^ x , y ). In the first 
node (m=S, m'=-S) H y °=H z =0 and H x °/T x =7l(n+l/2). If SH xy is small, it is easy to 
calculate y analytically: 

Y=A 2 [¥ x 2 +¥ y 2 +(7i 2 )(E DX 2 +E DY 2 )], (67) 
¥ x = 7i8H x /T x , ¥ y = 7iSH y /T y . (68) 
Thus, for E DX =E DY =E D (zero-field cooled after annealed at high temperature) exactly in 
the node Qi , x y =0) we get: 

Y-2ti 2 A 2 E d 2 . (69) 
When M in decrease, the width of the dipolar fields distribution increases. With 
increasing of E D the value of y in the nodes also increases. This result explains the 
experimental behavior of the tunneling splitting in the nodes (since As,-s is proportional 
to y m ). To check these formulas we have calculated y numerically using Eq. (12). The 
Figure 34 shows that (i) for different values of E D (i.e., depending on annealing) yin the 
nodes really behaves like E D 2 (see Figure (c)); (ii) around the nodes y has parabolic 
dependence on applied transverse magnetic field in accordance with Eq. (67) (see 
Figure (b)). The Figure 35 presents the same y as in Figure 34 (a) but calculated with a 
small "misalignment" angle m =l°. The curves in Figure 35 behave more similar to the 
experimental ones with respect to relative value of yin different nodes. 
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Above considerations concludes our discussion on the quantum regime of 
relaxation (low-T ground-state tunneling) in Feg and Mni 2 . Strictly speaking, all these 
effects are not restricted to these systems but are valid for all mesoscopic systems where 
quantum tunneling is associated with extremely small tunneling splitting, i.e. A~10" n 
with n much larger than one (n = 7 and n~ 11 in Fe 8 and Mni 2 ). The number of phonons 
available at these temperatures being negligible, phonons are not really relevant. On the 
contrary, other types of fluctuations are numerous in this regime and this is the case for 
fluctuations of the spin bath. The main result, which we would like to emphasize here, 
looks very simple and straightforward. All the physics in this limit depends on 
hyperfine and dipolar interactions. Any particular result is just the consequence of the 
first one. Finally we would like to say that the molecules Fe 8 and Mni 2 discussed in this 
paper are ferrimagnetic with a large non-compensated moment, i.e. their physics is 
dominated by the ferromagnetic order parameter. They have large spins and, therefore, 
important energy barriers and small tunneling splittings. In other systems with small 
non-compensated spins (eventually zero in antiferromagnetic molecules), energy 
barriers are much smaller leading to much larger tunnel splitting. In this case, the 
mesoscopic physics is not limited to the spin bath, the phonons bath is also very 
relevant. 
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